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ABSTRACT 

This contract supports the investigation of elemental abundances in the solar corona, 
principally through analysis of high-resolution soft X-ray spectra from the Flat Crystal 
Spectrometer on NASA’s Solar Maximum Mission. The goals of the study are a char- 
acterization of the mean values of relative abundances of elements accessible in the FCS 
data, and information on the extent and circumstances of their variability. This report is 
a summation of the data analysis and reporting activities which occurred since the last 
report, submitted two months early, in April 1994, to facilitate evaluation of the first 
year’s progress for contract renewal. Hence this report covers the period 15 April 1994 - 
15 December 1994. 


Subject terms: Solar corona, X-ray spectra, elemental abundances, abundance variability 



CORONAL ABUNDANCES AND THEIR VARIATION 


I. INTRODUCTION 

This is a semi-annual progress report for contract NASW-4814. The contract resulted 
from an award under NASA’s Supporting Research & Technology Program. Funding of 
a three-year contract began in June 1993. The contract was renewed for a second year 
in June 1994. The contract provides funding for the PI at about half-time effort. Since 
the last report was submitted two months early, in April 1994 rather than June 1994, to 
permit timely evaluation of progress and continuation of funding, this report covers an 
eight-month period, namely 15 April - 15 December 1994, rather than the nominal six- 
month period. Still, not as much was performed on this contract during this reporting 
period as originally planned, since the PI was heavily engaged in writing proposals over 
the summer and then in refereeing proposals (neither of which activity was charged to this 
contract). 

However, there were several significant accomplishments during the period, including a 
talk given at the joint American Geophysical Union/Solar Physics Division meeting in 
Baltimore in May, and an invited review at COSPAR meeting in Hamburg, for which a 
refereed paper is now in press in Advances in Space Research. Preparation for these talks 
and the feedback received at the meetings helped to clarify some ideas on several important 
issues involved in the contract research. Further, on the way to the COSPAR meeting in 
Hamburg, the PI visited the Service d’Atrophysique, Centre d’etudes de Saclay near Paris 
on her way to the COSPAR meeting in Hamburg and had some very useful discussions 
with Dr. Jean-Paul Meyer and Dr. Monique Arnaud on a variety of topics directly related 
to the research. (Saclay underwrote the costs for this portion of the foreign travel.) 

Also during this period, the PI has been working with colleagues Dr. Joan Schmelz (Rhodes 
College, Memphis, Tennessee) and Dr. Keith Strong (Lockheed, Palo Alto, California) on 
the draft of a paper which examines anomalous neon abundances in nonflaring active 
regions using a combination of two analysis techniques, as discussed below in §III.D.) 

II. SCIENTIFIC BACKGROUND 

The contract supports an investigation of elemental abundances in the outer atmosphere 
of the Sun, principally through analysis of high-resolution soft X-ray spectra from the Flat 
Crystal Spectrometer (FCS) on the Solar Maximum Mission ( SMM ), a NASA mission ded- 
icated to solar observations from 1980 through 1989. This instrument acquired an excellent 
data base for studying the relative amounts of oxygen, neon, magnesium, and iron in solar 
active regions in various states of evolution and activity. The project includes analysis of 
this data base to decouple the effects of temperature and abundance, to assess different 
theoretical calculations of spectral line intensities for use in the study, and to account for 
the possible effects of opacity due to resonance scattering of certain bright lines. The goals 
of the study are a characterization of the mean values of relative abundances of elements 
accessible in the FCS data, information on the extent of their variability, identification 
of possible correlations of variability with active region properties, and clarification of a 
possible association between abundance variability and active region dynamics. 
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Much of the observed variation in abundances between the photosphere and the corona 
appears to be associated with the first ionization potential (FIP) of the elements: elements 
with low FIP are systematically enhanced in the corona, solar wind, and solar energetic 
particles (SEP) relative to high-FIP elements, as compared with the photospheric com- 
position. Although there is no convincing detailed quantitative model yet available, the- 
oretical work in this area appears to be progressing, according to a review talk given by 
Jean- Claude Henoux at the COSPAR meeting in Hamburg. The basic idea is that some 
element separation mechanism operates in a temperature regime of about 10,000 K where 
low-FIP elements are ionized while high-FIP elements remain neutral. However, spectro- 
scopic data from SMM and other spaceraft show that coronal abundances in closed coronal 
structures present a more complex story than a simple two-step distribution of abundances 
(of high- vs. low-FIP elements) relative to photospheric values seen in the slow solar wind 
and in SEP events. There is also disagreement between the spectroscopic and particle 
measurements of the normalization of abundances relative to hydrogen, and in fact some 
disagreement among the different spectroscopic measurements. Knowing the abundances 
relative to hydrogen is crucial for many aspects of solar spectral data interpretation, be- 
ginning with converting line intensities to emission measure (i.e., the amount of emitting 
material). An attempt to reconcile some of the different spectroscopic measurements was 
made in the invited review of coronal spectroscopic results given by the PI in Hamburg in 
July 1994 (see §IV and attached preprint). 

Previous work by us (e.g., Strong et al. 1991, Saba & Strong 1993) and by others (e.g., 
McKenzie and Feldman 1992) has showed that ratios of low-FIP to high-FIP elements 
appear to vary at least about a factor of four, approximately from the photospheric values 
to the nominal coronal values defined by the SEP results. The precise upper and lower 
boundary values depend on the theoretical line emissivity calculations used for both the 
abundance line ratios and the temperature diagnostics. The question of the shortest ob- 
served timescale for variability hinges on the correctness of the line emissivity calculations 
as well, since the predicted emissivities are used to untangle changes in the abundance 
from evolution in the temperature. A key component in the predicted emissivity is the 
ionization fraction, i.e., how much of the element as a function of temperature is in the 
correct stage of ionization to emit the observed line; under typical “coronal conditions” of 
low density and local thermodynamic equilbrium, the ionization fraction is a factor which 
multiplies the effective probability of transition between the electron levels corresponding 
to the energy of the observed fine. 

III. CURRENT PROGRESS 

During the reporting period, 15 April to 15 December 1994, the PI has made progress in 
the following areas: 

(a) understanding better the issues and uncertainties in the ionization balance calculations 
used for the abundance analysis, 

(b) quantifying the role of resonance scattering in the measured line intensities, 

(c) examining measurements of absolute elemental abundances (i.e., abundances normal- 
ized relative to hydrogen), and 

(d) exploring anomalous neon abundances in active regions. 
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These are each discussed in more detail below. 

A. Ionization Balance 

The PI took advantage of her visit to Saclay to discuss at length with Dr. Monique Arnaud 
the uncertainties involved in ionization balance calculations in general, and in particular 
in the new calculations for iron by Arnaud & Raymond (I99x) [hereafter, ARA]. In some 
cases, the new calculations give significantly different results from previous calculations 
by Arnaud & Rothenflug (1985) [hereafter, ARO], and predict very different temperature 
behaviour for Fe XVII/Ne IX in the regime of 2-5 MK. The ARO calculations predict that 
this line ratio should be very flat over the temperature regime relevant for active regions, 
while the ARA calculations predict an increase with temperature, so that the deconvolution 
of derived abundance and temperature are substantially different. In comparisons with 
data from the SMM FCS, the new ARA calculations appear less consistent than the older 
ARO calculations for Fe XVII and Fe XVIII in the regime of a few MK. In particular, 
they imply that most of the time the relative Fe/Ne abundances in the FCS coronal 
measurements are less than the photospheric value, which would be very troubling. This 
would mean that the FlP-based picture is not even approximately right for spectroscopic 
measurements of closed coronal structured. However, it remains to be seen whether the 
apparent difficulties with the ARA calculations are in fact due instead to effects of resonant 
scattering or inadequate calculations for the Fe XVII or Fe XVIII excitation rates. It 
should be noted that, while the two calculations give different boundaries for the range 
°f Fe/Ne variability, the overall magnitude of the variability is comparable; also, while 
the ARA calculations imply that the variations on timescales less than an hour are due 
to temperature evolution rather than a change in abundance, day-to-day changes in the 
active region abundances remain. After discussions with Dr. Arnaud at Saclay, the PI has 
tentatively concluded that the best calculations would probably lie in between the ARA 
and ARO results for the cases relevant for the FCS analysis, so that these two extremes 
can be used to bound the uncertainties. The ionization balance calculations for may other 
abundant elements are currently being reconsidered by a community of theoreticians in 
preparation for SoHO analyses (H. Mason, priv. comm, at COSPAR meeting). 

B. Resonance Scattering Effects 

The recent paper by Waljeski et al. (1994) finds that the magnitude of the resonance 
scattering effect in active regions may be even larger for certain bright X-ray lines than 
had been found by Schmelz, Saba, and Strong (1992). Schmelz et al. suggested that the 
apparent depletion in the flux of the bright Fe XVII line at 15.01 A, a factor of two to 
three below predicted theoretical values compared to other bright resonance lines, might 
be due to resonance scattering. (This effect was discussed in detail in the last progress 
report.) The Fe XVII line at 15.0lA featured importantly in the early FCS abundance 
analysis. Waljeski et al. argue that it is insufficient to omit this one fine from the analysis, 
since the effect of resonant scattering on other lines used could also be significant. 

The new analysis by Waljeski et al. has the advantage of having considered two Fe XVII 
lines to estimate the importance of resonance scattering, the the 15.01 A line and a line at 
16.78A which has about 25% of the opacity of the 15.0lA line, and so it avoids possible 
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effects from uncertain relative abundances of different species. However, the estimated un- 
certainties of the predicted emissivities for the two Fe XVII line (A. Bhatia, priv. comm.) 
produce a large uncertainty in their expected line ratio, which in turn yields a large un- 
certainty in the derived opacity of the 15.0lA line, which is comparable to the difference 
between the two sets of results. When the PI analyzed additional data from the same 
active region, using another Fe XVII line which is essentially unaffected by resonance scat- 
tering, the magnitude of the effect was found to be closer to the original Schmelz et al 
result rather than the newer Waljeslci et al. result, even though in principal the latter was 
one more self-consistently than the former. That is, the optical depths for the Fe XVII 
line at 15 01 A appear to be closer to 2 than to 4, and the other lines used are essentially 
unaffected within the uncertainties. This interim result was reported at the AGU/SPD 
meeting in Baltimore in May, and will be incorporated in a paper being outlined which 
will compare the Fe XVII calculations with the FCS observationbs. 

Assessment of the magnitude of the effect of resonance scattering in the FCS data and 
the impact on the derived abundance variability continues. In the coming months, we 
will resume our comparison of the various theoretical calculations available for five bright 
Fe XVII lines covered by the FCS spectra and consult with several theoreticians who have 
performed some of the relavant calculations. Several theoreticians attending the COSPAR 
session on abundances (private comm, from A.H. Gabriel and F. Bely-Dubau) believe that 
the whole set of Fe XVII calculations should be revisited, to look for a basic problem with 
the calculation of the intensity of the lS.OlA line that is unrelated to resonance scattering 
or other transport effects. In the mean time, the abundance analysis is being continued 
using the Fe XVII line at 16.78 A which can be treated as optically thin to scatter. 

C. Absolute Elemental Abundances 

In the last few months, as a result of further thought and discussion, the PI has come to 
realize that, even if it is confirmed that resonance scattering is the sole source of the discrep- 
ancy between measured and predicted fluxes for the Fe XVII line at 15 . 0 lA, there are two 
problems with attempting to exploit resonance scattering to try to normalize abundances 
relative to hydrogen, as was attempted by Waljeski et al.: (i) the resonance scattering 
process in the active region geometry needs to be modeled more exactly to derive the 
calculated opacity to an accuracy better than about a factor of two; and (ii) the required 
substitution of the volume emission measure term for the integrated column density ne- 
glects a geometric term which is roughly proportional to L x f L /F A , where f L and F A are 
the volume and area filling factors, respectively. Thus the net systematic uncertainty could 
easily be comparable m magnitude to the typical low-FIP/high-FIP enhancement factor 
in the corona. These difficulties may be overcome with future observations which combine 
direct information on the appropriate geometric model (from high-resolution imaging) and 
independent information on the electron density, from a density-diagnostic line ratio which 

must be appropriate for the temperature regime of formation of the resonantly scattered 
line. 

Other methods of obtaining absolute abundances from spectroscopic measurements, and 
some attempts to reconcile the different results, were discussed in the COSPAR review 
paper and are noted in the accompanying preprint. In brief, in some cases it appears that 
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line-to-continuum flare abundance studies may have considered plasmas at temperatures 
sufficiently low that the free-bound component of the continuum was important, without 
using a self-consistent treatment of the elements contributing to the continuum. Some 
of the discrepancies between the SMM Bent Crystal Spectrometer and the Yohkoh Bragg 
Crystal Spectrometer abundance results may turn out to be related to differences in the two 
instruments. In particular, the Yohkoh BCS has a full-Sun field of view which may allow 
contamination of the “flare continuum” by cooler solar flux from regions not involved in 
the flare (J.L. Culhane, public comment during discussion at the COSPAR meeting). This 
would effectively reduce the derived abundances, which is consistent with the direction of 
the systematic SMM and Yohkoh discrepancy, although careful work by the Yohkoh BCS 
team members (many of whom were also involved in SMM abundance studies) is required 
to sort this out definitively. 

D. Neon Abundance Anomalies in Active Regions 

As discussed in the last report, two complementary methods of abundance analysis are now 
being combined to improve the FCS active region abundance study. The line-ratio analysis, 
which yields pair-wise relative abundances provided the emitting plasma is isothermal 
(or the emissivity functions have identical shapes in the relevant regime), is being used 
in parallel with a method which overlays curves of the line emissivities as a function of 
temperature. The different curves should overlap in a single temperature/emission measure 
solution if the plasma is indeed isothermal and the correct relative abundances are used 
or if reasonable adjustments to them can be made by simply shifting the curves vertically. 
An example of this second approach is shown in Figure 2 from the attached preprint of 
the COSPAR review paper. (The left panel shows nonintersecting curves for multithermal 
plasma during a flare; the right panel shows nearly isothermal plasma during the flare 
decay.) 

The advantages of combining this method with the current analysis are (i) that it provides 
a result that is less sensitive to any specific temperature diagnostic and (ii) that it makes 
it easier to keep track of the joint relative abundances of the four elements O, Ne, Mg, 
and Fe, without having to explicitly divide out the temperature dependences. In turn, this 
will guide the deconvolution of temperature and abundance effects which is the goal of the 
line-ratio analysis. 

As noted in the introduction section, the PI is now working with Drs. Schmelz and Strong 
to apply these combined methods to search for and confirm neon abundance anomalies in 
nonflaring active regions. A paper is being drafted which discusses the detection of anoma- 
lous neon abundances (with ratios of neon to oxygen both higher and lower than expected 
for photospheric values). Although the bulk of the FCS spectra show the expected ratios, 
within the statistic and systematic uncertainties, the occurrennce of anomalous abundances 
for the high-FIP neon, particularly in observations of quiescent active regions, are trou- 
bling for the commonly invoked FlP-based models of coronal composition. This anomalous 
variability shows that a competing element differentiation mechanism or mechanisms must 
apply. 

IV. PUBLICATIONS AND PRESENTATIONS 
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During the reporting period, the following presentations were made under this contract: 

A talk at the joint American Geophysical Union- American Astronomical Society/ 
Solar Physics Division Meeting in Baltimore, Maryland, in May 1994 was given on A 
Reexamination of the Impact of Resonance Scattering on FCS Active Region Abundance 
Measurements. 

An invited review on the subject Spectroscopic Measurements of Element Abundances 
in the Solar Corona: Variations on the FIP Theme was given at the 30th Assembly of the 
Committee on Space Research (COSPAR) held in Hamburg, Germany in July 1994, in a 
special one-day session on Element Abundance Variations in the Sun and the Heliosphere. 
A review paper summarizing the key points of the talk was submitted, refereed, and 
accepted for publication in Advances in Space Research, and a camera-ready manuscript 
was prepared. A preprint is attached. 

During this reporting period, two papers were published involving work done under the 
first year of the contract: 

The Composition of a Coronal Active Region, by K. Waljeski, D. Moses, K.P. Dere, 
J.L.R. Saba, K.T. Strong, D.F. Webb, and D.M. Zarro, appeared in the Astrophysical 
Journal in July 1994. A photocopy of the journal paper, to which the PI contributed 
substantially in the first year of the contract, is attached. 

The conference paper titled Implications of Coronal Abundance Variations, by Julia 
L.R. Saba and Keith T. Strong, appeared in July in the Proceedings of the Kofu Sym- 
posium, “New Look at the Sun with Emphasis on Advanced Observations of Coronal 
Dynamics and Flares - What Do We See with Yohkoh and Nobeyama Radioheliograph,” 
eds. S. Enome and T. Hirayama, Nobeyama Radio Observatory Report No. 360, 1994. A 
reprint is attached. 

A paper titled Anomalous Coronal Neon Abundances in Quiescent Active Regions is being 
prepared for submission to Astrophysical Journal Letters by J.T. Schmelz, J.L.R. Saba, 
D. Ghosh, and K.T. Strong. 

V. FOREIGN TRAVEL REPORT 

During the reporting period, Dr. Saba traveled to Hamburg, Germany, to attend the 
COSPAR (Committee on Space Research) meeting in Hamburg, Germany (11-21 July 
1994) and participate in two solar sessions: a one-day meeting (Session E2.1) on “Ele- 
ment Abundance Variations in the Sun and the Heliosphere,” on 16 July and a four-day 
symposium (Session E2.2) on “Solar Flare, Coronal, and Heliospheric Dynamics,” 18-21 
July. The first solar session, which Saba helped organize, brought together some of key 
researchers and the latest results in the rapidly changing field of solar abundances. The 
talks and the discussions following were extremely useful for understanding various as- 
pects of the coronal abundance picture and where some of the major uncertainties and 
discrepancies lie. 

At this session, Saba gave an invited review paper on spectroscopic measurements of coro- 
nal abundances including results from her FCS research to date. Questions and discussion 
after the talk suggested some areas which need more work or thought and some possible 
new approaches. Some of these have already been incorporated in the refereed conference 
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paper, which has been accepted for publication in Advances in Space Research; a preprint 
is attached. The other solar session, on dynamics, was pertinent to the next stage of 
Saba’s abundance research - an investigation of whether abundance variations are linked 
with coronal dynamics. This linkage is plausible and would mean that the observed abun- 
dance variations might be useful as a probe of conditions deeper in the solar atmosphere, 
but it remains to be demonstrated by the kind of correlative study planned. 

Saba was invited to visit the Service d’Astrophysique, Centre d’Etudes de Saclay, in Gif- 
sur-Yvette, near Paris, France during the week before the COSPAR meeting, with Saclay 
underwriting the expenses for the Paris leg of the trip. Saba was hosted there by Dr. 
Jean- Paul Meyer, who has done extensive work in the field of solar abundances. Saba and 
Meyer had many intensive, profitable discussions on topics such as approaches for integrat- 
ing the solar active region and flare abundance results, and reconciling the variable coronal 
abundances measured spectroscopically with the more stable coronal composition inferred 
from solar energetic particle measurements; the problem of normalizing abundances rel- 
ative to hydrogen; the roles played by elemental abundances in other plasma diagnostic 
measurements, and by line emissivity calculations in the inferred abundances. Some of the 
insights gained from these discussions were reflected in the COSPAR talk and subsequent 
paper. 

At Saclay, Saba also discussed with Dr. Monique Arnaud the issues involved in calculating 
ionization balances, and the problems in trying to use the new iron ionization fractions of 
Arnaud & Raymond (1985) in the FCS abundance analysis. The ionization fractions are an 
important factor in deriving temperatures and abundances from measured line intensities, 
and in deconvolving abundance effects from temperature evolution. 

Saba flew to Paris from Washington/Dulles on Thursday, 7 July 1994, arrived in Paris Fri- 
day morning, and worked at Saclay and at the home of Jean-Paul Meyer through Thursday, 
14 July. Thursday evening, she took an overnight train to Hamburg, arriving early Friday 
morning. She attended the COSPAR meeting through Thursday, 21 July, returned to 
Paris by overnight train Thursday evening, and then flew to Washington/Dulles on Friday, 
22 July. 
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ABSTRACT 

Solar wind and solar energetic particle (SEP) data yield systematic differences between elemental 
abundances in the corona and in the photosphere related to the first ionization potential (FIP) 
of the elements: low-FIP elements are preferentially enhanced relative to high-FIP elements by 
about a factor of four. Spectroscopic studies of the inner corona show that such a pattern may 
apply on average but not in detail for coronal loops: substantial abundance differences occur 
between different types of coronal structures, and variations have been found from flare to flare, 
from one active region to another, and over time in the same region; further, in some flares, 
anomalies such as enhanced Ne:0 ratios, distinctly at odds with the FIP pattern, show that a 
competing element selection mechanism sometimes operates. Details of the observed abundance 
variability - such as the magnitude of the variations, the relevant temporal and spatial scales, 
and correlations with other properties of the given coronal structure - may give important clues 
to the processes which supply and heat the corona, or they may reflect the changing physical 
conditions or locations where those processes take place. However, many such details remain to 
be established definitively. At present, abundance variability is primarily a major complication 
to data analysis and interpretation. However, once it is better understood, it may provide a new 
diagnostic tool for probing the lower layers of the solar atmosphere. 

INTRODUCTION 

Several recent comprehensive reviews of elemental abundances in the corona have been written 
/1 ,2,3,4/. The discussion here concentrates on several of the outstanding problems and incon- 
sistencies in spectroscopic determinations of coronal abundances. A detailed examination of the 
theoretical calculations and the atomic data needed to transform observed line intensities into 
abundances is given by H. Mason elsewhere in this issue. 

There has been fairly strong spectroscopic evidence for at least 20 years that coronal composition 
differs from that in the photosphere (e.g., the careful forbidden line study by Mason /5/). Yet 
many analyses of coronal abundances found results which appeared ambiguous or even contradic- 
tory. The discrepancies were due to a variety of sources: inadequacies in the atomic data, in the 
spectral or spatial resolution, or in incomplete separation of temperature and abundance effects; 
in some cases, the uncertainties were simply too large for any definitive statements to be made. 
One possibility that was not at first considered was abundance variability within the corona, 
discovered for flares by Sylwester, Lemen, and Mewe /6/, and then for active regions by Strong 
and colleagues /7,8/; variations in composition from one type of coronal structure to another 
were discovered by Widing and Feldman /9/ and colleagues (see /4 / and references therein). 

It is currently difficult to synthesize the existing spectroscopic results into a single coherent picture 
when some of the most intriguing results appear contradictory and important questions remain 
to be answered. Still, some broad statements may be made: 

(1) There appears to be basic agreement with the SEP and solar wind picture that elemental com- 
position in the corona differs on average from that in the photosphere, with an overall systematic 
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bias based on first ionization potential (FIP) of the elements, with low-FIP (<10 eV) elements 
enhanced relative to high-FIP (>11 eV) elements, compared to their photospheric ratios. 

(2) The case for abundance variability within the corona is very strong. This variability includes 
systematic variations in average composition from one type of coronal feature to another /9/, 
variations between flares (e.g., /6,10,11/), variations between active regions /7, 8, 12,13/, and 
variability over time in a given active region /8,13/. 

In many cases, the variations are consistent with some generalized version of the FlP-based 
picture, if one allows for variation in the boundary conditions for the element separation mecha- 
nism^), or perhaps for different heights from which source material is drawn into the corona. In 
some cases, the abundance variability includes effects not consistent with the general FlP-related 
bias; notably this includes enhancements of neon (a high-FIP eiement) during certain flares and 
long-duration events /14,15,16/, which Shemi / 1 7/ explains as being due to preferential pho- 
toionization of neon. There is growing spectroscopic evidence that the FIP enhancement factor, 
rather than being a simple, constant step function, may be a more complicated function /4,18/, 
or that it may be variable /7, 12, 13, 19, 20/. Recent analyses of active region abundances using 
soft X-ray spectrometer observations from P78-1 /12/ and SMM /13/ suggest that the “coronal” 
ratios of low-FIP/high-FIP element pairs such as Fe/Ne, Fe/O, Mg/Ne, and Mg/O can vary 
approximately from their photospheric values to the “SEP coronal” values or even higher. 

ABUNDANCE ANALYSIS METHODS 

Several kinds of spectroscopic techniques are used to determine coronal abundances. The most 
commonly used methods and some of their difficulties are discussed here briefly. 



Fe XVIII / Fe XVII 

Fig. 1. Abundance-diagnostic ratio Fe XVII/Ne IX plottec vs. temperature parameter, 

Fe XVIII/Fe XVII. (From Saba and Strong /13/.) 

Line Ratios 

One technique for determining relative elemental abundances is to use abundance-independent, 
temperature-diagnostic line ratios, such as the ratios of lines from different ionization stages of 
iron, to characterize the temperature of the plasma, and then to examine in parallel the flux ratios 
of lines from different elements, which are temperature-invariant or which vary only slowly with 
temperature over the relevant temperature regime. An example of this method’s use in analysis 
of FCS data from SMM / 13/ is shown in Figure 1. Different symbols designate different active 
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regions observed. Two sequences of measurements taken in the decay phase of long-duration 
events are connected by solid curves. The line ratios can be converted to abundance ratios by 
dividing out the temperature dependences of the line emissivity functions. The two different 
theoretical emissivity curves plotted (dashed and dotted curves), discussed below, imply different 
conversion algorithms. This method was also used by McKenzie and Feldman / 12/, for data from 
P78-1, and in a series of papers by Widing and Feldman /9, 19,21/ and colleagues, for analysis of 
EUV data from Skylab. 

The line-ratio technique is dangerous to use blindly, if the shapes of the emissivity functions of 
the two lines being used for the abundance-diagnostic ratio do not match with good precision 
in the temperature regime where there is significant emission measure; however, this method is 
convenient to use and complementary to other methods. In particular, line ratios can provide esti- 
mates of relative abundances for use in emissivity curve overlay and for a first guess in differential 
emission measure (DEM) analyses; consistency checks on the line-ratio method can be made by 
comparing the results from multiple line ratios, for lines covering a broad range of temperature. 

Emissivity Curve Overlay 

A more general method than the line-ratio technique is to overlay the emissivity functions for a 
variety of observed emission lines. This method was applied to soft X-ray data by Strong /22/, and 
further developed by Veck et al. /23/. If the isothermal approximation is valid, then the curves 
will intersect in a single emission-measure/temperature solution, to within the uncertainties, 
provided that the theoretical emissivities are correct and a correct set of abundances is adopted. 
If the plasma is not isothermal to within the measurement uncertainties, or if an incorrect set 
of abundances is assumed, the curves will not all overlap in a single location. Figure 2, taken 
from Veck et al. /23/, shows examples of these two cases. For serious quantitative studies, this 
method should be accompanied by a rigorous analysis of the uncertainties. 




Fig. 2. Overlay of emissivity curves from soft X-ray lines. The left panel shows curves 
for multithermal plasma during a flare; the right panel shows nearly isothermal plasma 
during the flare decay. (From Veck et al. /23/.) 


The emissivity curve overlay method has been used recently to look for abundance anomalies in 
SMM FCS measurements of active regions, for a set of active regions (see paper by J. Schmelz 
in this issue), and to determine the composition of a single active region for which soft X-ray 
broadband data were available in addition to soft X-ray emission line data /24 / ; the emissivity 
curves for the broadband data were derived using synthetic spectra at various temperatures for 
a given set of abundances. In the latter study, Waljeski et al. /24/, compared the results for 
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isothermal models and for standard empirical Gaussian differential emission measure (DEM) 
models; the abundance results were found to be indistinguishable. 

Differential Emission Measure 

If there are enough emission lines, with good temperature coverage and sufficient overlap, it is 
possible to do a DEM analysis of temperature, which involves an iterative inversion procedure 
/25/. Since the elemental abundances are not known a priori , it is useful to have several lines 
from different ionization stages of the same element. If the emission lines from other elements can 
be scaled to match the DEM at the appropriate temperature, the relative elemental abundances 
can be determined from the normalization adjustment needed to make the different lines mutually 
consistent. An incorrect set of abundances, or continuing the iteration procedure past the point 
where the model falls within the one-sigma uncertainties of the data, can introduce artificial 
complexity into the DEM solution / 1 8/. In general, a full DEM treatment for many spectra 
may be impractical, but it is possible to examine several cases for general guidance in the data 
analysis, and to compare the results with synthetic multithermal spectra calculated from first 
principles for various reasonable physical models. 

Line-to-Continuum Studies 

This method has been applied almost exclusively to soft X-ray flare studies. If used correctly, 
for a detector with a low intrinsic background, line-to-continuum intensity ratios can be used to 
obtain elemental abundances relative to hydrogen. However, there are several additional caveats: 
The method requires that the emission line and the adjacent (contamination-free) continuum 
be formed in the same volume; that the plasma be in ionization equilibrium, with a Maxwellian 
electron velocity distribution; and that the free-bound contribution to the continuum (which itself 
depends on the coronal composition) be accounted for properly. This last requirement is difficult 
to meet at temperatures below about 15 MK for wavelengths shorter than about 10 A, without 
a priori knowledge of the abundances of certain elements (see, e.g., study by Culhane /26 /). 
On the other hand, the method is most accurate for plasmas which are at least approximately 
isothermal, a condition not typical of flares. Both the free-free and the free-bound continuua 
depend significantly on the abundance of helium, w*hich is generally unknown or known poorly in 
the corona (however, see paper by J.L. Culhane et al. in this issue). 

The continuum intensity (photon s^A"* 1 ) emitted by an isothermal (Maxwellian) plasma is 
approximately given by: 

Ic(T) = JiGjj + Gjt + G^T-^X-^-^n^dV, (1) 

where the average Gaunt factors G//, G/j, G 27 represent contributions to the continuum emission 
from free- free, free-bound, and two-photon processes. The continuum Gaunt factors and the 
resulting continuum have been calculated by various authors; the calculations by Mewe, Lemen, 
and van den Oord /27/, which use the cosmic abundances of Allen /28/, are given in a convenient 
form and have been used in most of the SMM and Yohkoh analyses. Figure 3 shows plots from 
Mewe et al. giving the continuum components at plasma temperatures of 10 MK and 20 MK. 

At 20 MK, the free-free component dominates; at active region temperatures, the free-bound 
component is important for wavelengths as large as about 20 A /26/. At temperatures near 10 
MK, for Allen abundances, the free- free and free- bound components are comparable at a few 
A /27/. Hence the line-to-continuum method must be used witii caution, with a self-consistent 
treatment of the abundances in the free-bound component, for flare temperatures below about 
15 MK and wavelengths below about 10 A, and perhaps even at higher temperatures and longer 
wavelengths. If the elements contributing to the local free-bound emission are enhanced (relative 
to hydrogen) compared to cosmic composition, the free-bound component will be higher than 
anticipated, leading to artificially low line-to-continuum intensity ratios and hence artificially low 
abundances relative to hydrogen. 
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In concentrating on flare decays, in an effort to avoid transient nonequilibrium conditions and 
minimize nonthermal motions, some SMM and Yohkoh line- to- continuum studies (e.g., /29/) have 
included plasma at temperatures where the free-bound component is important, and so may have 
introduced significant systematic errors. The same problem might also have affected some of the 
results of Veck and Parkinson /30/. 


10 * 10 * 10 * 10 ® 10 * 10 * 



Fig. 3. Relative contributions of free-free (dotted curve), free-bound (dot-dash curve), 
and two-photon (dashed curve) to the total continuum emission (solid curve) at 10 MK 
(left) and 20 MK (right), for Allen abundances. (From Mewe et al. /27/.) 

There appears to be a systematic difference between the calcium abundance results from SMM 
and from Yohkoh , with the Yohkoh values falling preferentially at the low end of the SMM values, 
and showing much lower overall variability. Potential sources of this difference are being explored 
by the Yohkoh BCS instrument team members, many of whom were also previously involved in 
SMM BCS analyses. One possibility which is being explored is that the “flare continuum” in the 
Yohkoh BCS instrument is contaminated by cooler solar flux from regions not involved in the flare 
/31 /; one of the main differences between the SMM and Yohkoh BCS instruments is the differing 
fields of view - the SMM BCS has a 6 arc minute (FWHM) field of view, while the Yohkoh BCS 
views the whole Sun. 

Optical Forbidden-Line Equivalent Width 

Pottasch /32/ describes the procedure for determining absolute abundances from the equivalent 
width of optical coronal forbidden lines. Conceptually, this method is similar to the line-to- 
continuum method described above, but the implementation and caveats are different. The 
visible continuum in the lower corona is produced almost entirely by electron scattering of pho- 
tospheric radiation /32/. Hence the primary concerns are determining the electron density in the 
same region as that where the equivalent width is measured, and removing the effects of line-of- 
sight averaging of the emission line in the optically thin corona. One uses statistical equilibrium 
equations to calculate the number of ions in all levels from the population of the upper level, 
which requires accurate transition probabilities and electron collisional excitation rates. Then 
one sums all of the different ionization stages of the given element, which requires knowledge of 
the temperature distribution and accurate ionization balance calculations. Using this method, 
Pottasch /32 / found an abundance ratio of Fe/H = 8 x!0~ 5 , about 2.5 times the current estimate 
of the photospheric value. With greatly improved atomic data from detailed calculations includ- 
ing cascades, Mason /5/ found similar coronal enhancements of Ca/H and Fe/H in a coronal 
condensation. 
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Generaliz e d DEM 


Fludra et al. /1 1 / proposed a method which combines the well known technique of using line 
ratios of the Ca XIX and Fe XXV resonance lines w and their dielectronic satellites k and j , 
respectively /33/, with standard multi-temperature analysis /34/; the generalized method uses 
line fluxes and line flux ratios simultaneously to derive the DEM to obtain relative abundances or 
the absolute abundance if continuum flux for one of the lines is also included in the analysis. This 
method is especially useful for obtaining the Fe/Ca and the Fe/H abundance in flares observed 
by the SMM BCS. (A background resulting from fluorescence of germanium crystals contributes 
a significant, time- variable proportion of the total background in the SMM BCS iron channel, 
while the SMM Ca XIX background is generally regarded as free of contamination.) 

ASSUMPTIONS 

A variety of assumptions go into abundance analyses of coronal spectroscopic data. Generally the 
plasma is assumed to be optically thin, at least to absorption, and typically to electron scattering 
as well. The plasma is assumed to be in ionization equilibrium; in flare studies this assumption 
may not be valid in the impulsive phase, which is therefore usually excluded from the analysis. 
The electron velocity distribution is assumed to be Maxwellian. Whenever a line ratio or line-to- 
continuum ratio is taken, an implicit assumption is made that the emissions being compared come 
from the same volume(s). On the other hand, Jakimiec et al. /35/ point out that the volume of 
plasma which produces flare continuum emission may be larger than that which produces flare 
line emission if the electron mean free path is large enough to allow electrons to escape from the 
high-temperature volume before interacting with the highly ionized calcium and iron atoms. In 
effect, this means that the electron distribution which produces the line emission might have a 
depleted high-energy tail, which would result in lower observed values of coronal abundances. 

Theoretical emissivitv calculations 

To convert relative line intensities to relative abundances or line-to-continuum or equivalent-width 
measurements to absolute abundances requires modeling of the emission line and continuum 
emissivities as a function of temperature and knowledge of the temperature distribution of the 
emitting plasma. The theoretical calculations of emissivities are discussed in detail in a review 
elsewhere in this issue. Here some general aspects of the calculations are discussed. 

In the “coronal approximation” of a low-density plasma in thermodynamic equilibrium, the emis- 
sivity function of an emission line is essentially a product of the elemental abundance times two 
temperature-dependent factors: the fraction of the emitting element in the relevant ionization 
stage (the ionization fraction) and the probability that the ion will make the given transition 
to produce the line (which depends on the level population and the excitation rate). For the 
H-like and He-like lines, the level populations and excitation rates are believed to be well under- 
stood. Multielectron ions, such as the neon-like Fe XVII and fluorine-like Fe XVIII resonance 
lines are harder to model, but the excitation rates are believed to be knowrn to about 30% /36/. 
The ionization fractions, especially for multi-electron ions, are more difficult to calculate and 
correspondingly more uncertain. 

A variety of calculations of equilibrium ionization balance have been made; differences in the 
calculations change the slope and the absolute scaling of the net line emissivity functions, which 
in turn affect the theoretical line or line-to-continuum ratios and hence the derived values of 
abundances. If the calculated ionization balance is in error, the derived abundances may show 
some apparent temperature dependence. A lack of ionization equilibrium would also change the 
ionization fractions. However, currently accepted flare models and densities point to ionization 
equilibrium during most of a flare, except perhaps the early rise phase /37/. Disagreement of 
various existing ionization balance calculations for calcium or iron with observed X-ray spectra 
has been reported previously by many authors, and empirical corrections have been proposed by 
some of them (e.g., /10, 1 1 ,38,39/) . None of the empirical corrections explain fully all observed 
line-to-line or line-to-continuum ratios, and they should be regarded with caution. 
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Recent calculations of iron ionization fractions by Arnaud and Raymond /40/ [hereafter, ARA] in 
some cases give significantly different results from previous calculations by Arnaud and Rothenflug 
/41/ [hereafter, ARO], as shown by the two theoretical curves in Figure 1, which predict very 
different temperature behaviour for Fe XVII/Ne IX. The ARO calculations predict that this line 
ratio should be very flat over the temperature regime relevant for active regions, while the ARA 
calculations predict an increase with temperature. In comparisons with data from the SMM FCS, 
the new ARA calculations appear less consistent than the older ARO calculations for Fe XVII and 
Fe XVIII in the regime of a few MK / 13/. However, it remains to be seen whether the apparent 
difficulties with the ARA calculations are in fact due instead to effects of resonant scattering 
or inadequate calculations for the Fe XVII or Fe XVIII excitation rates. It should be noted 
that, while the two calculations give different boundaries for the range of Fe/Ne variability, the 
overall magnitude of the variability is comparable; also, while the ARA calculations imply that 
the variations on timescales less than an hour are due to temperature evolution rather than a 
change in abundance, day-to-day changes in the active region abundances remain. 

Optically thin emission 

Although coronal emission lines are generally thought to be optically thin, results from a recent 
reconsideration of radiative transfer of Fe XVII in the corona suggest that resonance scattering 
might be responsible for the observed depletion of flux (relative to other bright emission lines) in 
the Fe XVII resonance line at 15.01 A in active regions /42/. 

In resonance scattering, a resonance line photon is absorbed by an ion in the ground state and then 
re-emitted, generally in a different direction. This absorption and re-emission is indistinguishable 
from scattering and depends on the geometry of the region being observed. In general, photons 
would be preferentially scattered out of the line of sight for relatively dense, bright areas, as 
compared to fainter surrounding areas. The absorption coefficient for resonance scattering, as 
given by Acton /43/ for J = 1 — ► 0 transitions, is 

k 0 = 9.31 x 10 - lB fi u AiAn t \ A M lf2 T? l/3 , (2) 

where // u is the absorption oscillator strength, A{ is the temperature-dependent relative abun- 
dance of the ion, A is the elemental abundance relative to hydrogen, n e is the density, is the 
wavelength in A, M is the mass in atomic mass units, and T t is the ion temperature in units of 
10 6 K. The corresponding opacity at line center is 

T’o = J KodL, (3) 

where dL is the line-of-sight element. If resonance scattering is in fact responsible for the apparent 
flux reduction of the Fe XVII line at 15.01 A, then it must be taken into account in deriving relative 
abundances from the measured line flux ratios, or other unaffected lines must be used to avoid 
systematic effects in the derived values of relative abundances. It is important to determine if 
resonance scattering does play a major role in the FCS and in other data sets for reasonable 
coronal parameters. 

NORMALIZATION OF ABUNDANCES 

A key issue is whether the low-FIP elements are enhanced or the high-FIP elements are depleted 
(or both) in the corona relative to their photospheric values. The absolute normalization of abun- 
dances is important both for determining the amount of source material producing a given line 
intensity and also for modeling the element separation mechanism(s) responsible for the differ- 
ences between photospheric and coronal composition. Since the bulk of the plasma is hydrogen, 
the question hinges on the behavior of elements relative to hydrogen. The solar wind and SEP 
observations seem to imply that the low-FIP elements are enhanced in the corona by a factor of 
3 to 4 while the high-FIP elements have their photospheric values. A variety of acceleration and 
transport processes may treat protons (i.e., ionized hydrogen) differently than heavier ions, so 
that heavy ion-to-proton ratios may not directly reflect the coronal heavy element-t.o-hydrogen 
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composition. Hence, it is desirable to have independent information from the spectroscopic abun- 
dance measurements. Unfortunately, at present, the spectroscopic results are not consistent. 

In the upper chromosphere, transition region, and corona, the temperature is above 3 x 10 4 
K, hence hydrogen is essentially fully ionized, and therefore has no emission lines. Thus, most 
abundance studies of emission lines from these regions give their results in terms of relative 
abundances. Recent EUV spectroscopic abundance results have been interpreted to mean that 
the low-FIP elements are enhanced exactly as in the SEP/soiar wind picture /4/. However, 
while the arguments are plausible, they are somewhat indirect and involve comparison of lines 
from different temperature regimes. Line-to-continuum and equivalent-width analyses give direct 
measures of abundances relative to hydrogen, subject to the caveats listed above. However, the 
line-to-continuum measurements show a large range of results, from a factor of 4 enhancement 
of low-FIP elements / 44/, to factor of two or less enhancement of low-FIPs, accompanied by a 
factor of two depletion of high-FIPs / 18,30/. Some of these measurements may be compromised 
by inadequate treatment of the free-bound continuum, or by contamination of the continuum 
by solar or nonsolar flux. Equivalent-width studies have shown factor of 2 to 3 enhancement of 
low-FIP elements /5,32/. Now that there are better atomic data available, it would be useful to 
examine more of the eclipse forbidden line data sets. 

Two novel approaches have been suggested recently for determining absolute coronal abundances. 
In the first of these (see paper by K. Phillips et al. elsewhere in this issue), Phillips et al. /45,46/ 
compared the intensities of Fe Ka and K (3 lines (which they argue are formed principally by 
fluorescence in the photosphere) to the Fe XXV resonance line (produced by hot flare plasma 
in the low corona). They deduced that the Fe/H coronal abundance ratio is not more than a 
factor of two higher than photospheric for a large sample of different kinds of flares observed by 
SMM BCS, P78-1 SOLEX or Yohkoh BCS. It is not clear if this result could be affected by an 
underestimation of the free-bound component of the continuum emission, as discussed above for 
line-to-continuum studies. The possibility should also be considered that the Fe Ka and K/3 lines 
are formed at heights considerably above the photosphere, where the element differentiation has 
already occurred. Still, if these potential complications can be ruled out, this result implies that 
the low-FIP elements are not greatly enhanced in the corona. 

The second new* method proposes to exploit resonance scattering of a bright Fe XVII line at 15.01 
A to get absolute abundances in active regions. With this method, using SMM FCS emission line 
intensities and AS&E broadband X-ray image data for information on geometry, Waljeski et al. 
/24/ found Fe:H = 3.2 x 10~ 4 (±60%), about 10 times the phoiospheric value. In this method, 
the effective optical depth at line center, r, to scattering is determined from the observed flux 
depletions of a pair of the Fe XVII lines and then also calculated from equations (2) and (3). For 
a line unaffected by resonance scattering, the line intensity is directly proportional to the factor 
Fe/H times the volume emission measure, EM = / n t 2 dV , where dV is the elemental volume 
(times factors involving atomic data and constants). For unifDrm density, / n e dl ~ n c L and 
/ n^dV n e 2 La, w*here a is the cross-sectional area and L is the path length* which can be 
estimated from imaging data. Making the substitution / n t dL zt [f n*dV]/[n e a] in the equation 
for r allows one to solve for Fe/H. 

Unfortunately, even if it is confirmed that resonance scattering is the sole source of the discrep- 
ancy between measured and predicted fluxes for the Fe XVII line at 15.01 A, there are currently 
two problems with using this approach: (1) the resonance scattering process in the active region 
geometry needs to be modeled more exactly to derive the calculated opacity to an accuracy better 
than about a factor of two; and (2) the required substitution of the volume emission measure 
term for the integrated column density neglects a geometric term which is roughly proportional 
to L x /i//>t, where fi and f A are the volume and area filling factors, respectively. Thus the 
net systematic uncertainty could easily be comparable in magnitude to the low-FIP/high-FIP 
enhancement factor in the corona. These difficulties may be be overcome with future observa- 
tions which combine direct information on the appropriate geometric model (from high-resolution 
imaging) and independent information on the electron density, from a density-diagnostic line ratio 
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(which must be appropriate for the temperature regime of formation of the resonantly scattered 
line). 

CONCLUSIONS 

Although a large amount of recent work has led to many interesting new results, the state of 
spectroscopic measurements of coronal abundances is not as mature as that of the solar energetic 
particle (SEP) and solar wind measurements. There remain ambiguities and contradictions in the 
spectroscopic coronal abundance results. The basic uncertainties include: discrepancies between 
different instruments and different analysis techiques on absolute abundances, and discrepan- 
cies between different instruments on the degree of abundance variability between and within 
flares. Our knowledge is currently sparse in several important areas: the relevant temporal scales 
and spatial scales for abundance variability in nonflaring regions, and the degree and nature of 
variability in low'-FIP/low-FIP and high-FIP/high-FIP pairs of elements. We still lack basic ob- 
servational information on the distribution of abundances with height. We do not yet know how 
to relate the spectroscopic and SEP results. All of this information is needed to provide realistic 
constraints on models of the element selection mechanism(s) responsible for the differences in 
composition between the photosphere and the corona. 

Although we are still disentangling the detailed effects of variable abundances from the effects 
of uncertainties in the emissivity calculations and physical modeling, it is clear that significant 
coronal abundance variations occur. There are many cases for which the apparent variations 
are much larger than the uncertainties, which cannot be explained in terms of problems with 
emissivity calculations or the presence of multithermal plasma. Continuing work on the theoretical 
predictions for line emissivities (particularly the ionization fractions), on accurate modeling of 
continuum emission (particularly the abundance dependence of the free-bound component), and 
on reasonable physical modeling of the geometry and multithermal distribution of plasma is 
needed to derive the details of the abundance variations required to challenge physical models. 
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ABSTRACT 

The relative abundances of iron, oxygen, magnesium, and neon in a coronal active region are determined 
from measurements of soft X-ray line and broadband intensities. The emission measure, temperature, and 
column density are derived from these measured intensities and are used to place a constraint on the abun- 
dances of the heavier elements relative to hydrogen in the corona. 

The intensity measurements were made on 1987 December 11, when an active region was observed jointly 
by the American Science and Engineering (AS&E) High Resolution Soft X-Ray Imaging Sounding-Rocket 
Payload and the X-Ray Polychromator Flat Crystal Spectrometer (FCS) onboard the Solar Maximum Mission 
spacecraft. The coordinated observations include images through two broadband filters (8-29 A and 8-39, 
44-60 A) and profiles of six emission lines: Fe xvn (15.01 A), Fe xvu (15.26 A), O vm (18.97 A), Mg xi 
(9.17 A), Ne ix (13.44 A), and Fe xvin (14.21 A). 

The effects of resonance scattering are considered in the interpretation of the FCS line intensities. We calcu- 
lated the expected intensity ratio of the two Fe xvii lines as a function of optical depth and compared this 
ratio with the observed intensity ratio to obtain the optical depths of each of the Ifnes and the column density. 

The line intensities and the broadband filtered images are consistent with the emission from a thermal 
plasma where Fe, O, Mg, and Ne have the “adopted coronal" abundances of Meyer (1985b) relative to one 
another, but are not consistent with the emission from a plasma having photospheric abundances: The ratios 
of the abundances of the low first ionization potential (FIP) elements (Fe and Mg) to the abundances of the 
high-FIP elements (Ne and O) are higher than the ratios seen in the photosphere by a factor of about 3.5. 
This conclusion is independent of the assumption of either an isothermal or a multithermal plasma. 

The column densities derived from the Fe xvn line ratio and the geometry of the active region provide a 
diagnostic of the abundance of hydrogen relative to the heavier elements. We find that the abundance of iron 
with respect to hydrogen in the corona is higher than the value given Meyer (1985b) by a factor of 8.2 (-5.4, 
+ 5.1). This means that, for the observed active region, the absolute abundances of the low-FIP elements (Fe 
and Mg) are enhanced in the corona relative to the photosphere, while the abundances of the high-FIP ele- 
ments (Ne and O) are either slightly enhanced in the corona or nearly the same in the photosphere and the 
corona. 

Subject headings: Sun: abundances — Sun: activity — Sun: X-rays, gamma rays 


1. INTRODUCTION 

Knowledge of the composition of the corona is an important 
key to our understanding of the physical processes of the solar 
atmosphere. The coronal elemental abundances reflect in some 
manner the way in which solar plasma is transported, acceler- 
ated, and heated. The composition of the corona must also be 
known in order to interpret observations of coronal emission: 
The relative abundances of the heavy elements must be know n 
to determine temperature from the intensity ratios of coronal 
lines of different elements, w hile the absolute abundances (the 
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abundances of the heavier elements with respect to hydrogen) 
must be known to determine emission measure and radiation 
loss. 

In this study, broadband images and emission intensities of a 
coronal active region are used to constrain the relative and 
absolute abundances of Fe, O, Mg, and Ne. The relative abun- 
dances are measured by comparing the observed line and 
broadband intensities with the intensities that would be 
expected from model plasmas of a specific composition. The 
optical depths of the Fe xvn lines and the geometry of the 
region provide the additional information necessary to deter- 
mine absolute abundances. 

1.1. Relative Abundances of Fe. Mg. O. and Ne 
There is considerable evidence, from in situ solar wind 
(reviewed b> Geiss 1982) and solar energetic particle measure- 
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ments (e.g.. Brennernan & Slone ] 9S5 : Meyer 1985a). as well as 
from coronal spccircscopic measurements (e.g.. Yeck & 
Parkinson 1981: Wiaing & Feldman 1989). that ihe composi- 
tions of the solar corona and the solar wind differ in a system- 
atic way from the composition of the photosphere. Meyer 
i l^sfb) synthesized the abundance measurements published by 
a large number of investigators. On the basis of this observa- 
tional evidence. Meyer (1985b) proposed that the composition 
of the corona and solar wind is fractionated according to first 
ionization potential (F1P), so that the ratio of the abundances 
of the low-FIP elements to the abundances of the high- FI P 
elements is higher in the corona and the solar wind than in the 
photosphere. Meyer (1985b) presented a set of “adopted 
coronal" abundances, which differed from the photospheric 
abundances in that the ratios of the abundances of the low-FIP 
elements to the abundances of the high-FlP elements were 
higher in the corona than in the photosphere by a factor of 
about 3.5. 

In this paper, we compare the intensity observations to the 
emission expected from plasmas with two sets of relative abun- 
dances: photospheric abundances (from Grevesse & Anders 
1991 with the more recently determined photospheric iron 
abundance from Grevesse, Noels, & Sauval 1992) and F1P- 
fractionated “adopted coronal” abundances (Meyer 1985b). 
Although the observed active-region plasma is contained by 
closed magnetic loops that are connected to the photosphere, 
we find that the observed intensities are not consistent with a 
plasma having photospheric abundances. The observations are 
consistent with thermal plasmas having the FlP-fractionated 
adopted coronal relative abundances of Meyer (1985b). 

Previous studies of Flat Crystal Spectrometer (FCS) data 
(e.g., Nitta et al. 1991; Saba & Strong 1991) have shown that 
some active-region spectra were more consistent with Meyer's 
adopted coronal relative abundances than with photospheric 
abundances; however, the present comparison of FCS spectra 
and broadband filter data gives a more rigorous quantitative 
basis for this conclusion. 

1 .2. A bsolut e A bundances 

Once the relative abundances of the high-FIP and low-FIP 
heavy elements have been established, a measurement of the 
hydrogen abundance is necessary to determine whether the 
high-FIP elements are underabundant in the corona relative to 
the photosphere or whether the low-FIP elements are over- 
abundant. On the basis of spectroscopic observations of line 
and continuum emission from the corona (Veck & Parkinson 
1981), Meyer (1985b) concluded that the abundances of the 
high-FIP elements were depleted in the corona relative to the 
photosphere. However, recent evidence, reviewed by Meyer 
1993a, b) and Feldman (1992), suggests that‘the abundance of 
hydrogen relative to the heavier elements is lower than predict- 
ed by Meyer (1985b), so that the low-FIP elements are over- 
abundant in the corona. 

There is no separate measurement of the continuum in our 
set of line and broadband observations from which to directly 
measure the abundance of hydrogen. The broadband filler- 
grams have the highest contribution of continuum to the 
signal. However, for all of the filtergram images, the continuum 
is expected to contribute less than 20% of the total filtered 
spectrum for emission from plasma at active-region tem- 
peratures with the adopted coronal abundances (Mever 
1985b). 

The effects of optical depth on the Fe xvn lines are used here 


as a diagnostic of absolute abundance. Most previous studies 
have assumed that the corona is optically thin to soft X-rays. 
As Acton (1978), Rugce & McKenzie (1985), and Schmelz, 
Saba. & Strong (1992) suggested, we find that the optical depth 
of the lines must be considered in the analysis of these broad- 
band and line intensities. The intensity ratio of the two Fc xvn 
lines (15.01 A/15.26 A), which differ significantly in oscillator 
stength, is used as an optical depth diagnostic. Significant 
optical depths are found for these two Fe xvn lines, the O vm 
line (18.97 A), and the Ne tx line (13.44 A). The optical depth of 
the Fe xvn line at 15.01 A has been proposed as a density 
diagnostic by Schmelz et al. (1992). Here, we use the column 
densities, J n t dl , together with the emission measures, J n] dl , of 
coronal soft X-ray lines to provide information on the abun- 
dance of hydrogen. 

The abundances of the heavier elements relative to hydrogen 
are found to be a factor of 8.2 ( + 5.4, -5.0) higher than the 
abundances tabulated by Meyer (1985b), implying that hydro- 
gen behaves as a high-FIP element in the corona. Therefore, 
the abundances of the low-FIP elements are enhanced in the 
corona relative to the photosphere, while the abundances of 
the high-FIP elements, including hydrogen, are nearly the 
same. 

2. OBSERVATIONS 

On 1987 December 13, a solar active region, with NOAA 
designation AR 4901, was observed simultaneously by the 
American Science and Engineering (AS&E) High Resolution 
Imaging Sounding-Rocket Payload and the X-Ray Poly- 
chromator Flat Crystal Spectrometer (FCS) onboard the Solar 
Maximum Mission ( SMM ) spacecraft. The active region was 
located near the southwest limb of the Sun, at approximated 
S22 W72. 

2.1. Spectrometer Observations 

The FCS, described by Acton et al. (1980), consisted of seven 
rotatable Bragg crystal spectrometers that were able to simul- 
taneously observe the same field of view. The 1987 December 
11 observations were all obtained with the first and third 
crystal channels. 

The response of the FCS collimator to an infinite uniform- 
brightness source was approximately pyramidal, with a full 
width of 30*. Each flux measurement made by the FCS is the 
convolution of this field response function with the intensity of 
the source. 

Two types of observations were made during the SMM 
satellite orbit (17:55 UT to 38:30 UT) of the AS&E sounding- 
rocket flight: maps of line peak intensity and spectral scans of 
line profiles. The maps were made by holding the angle of 
incidence of the FCS crystals fixed, to select a wavelength at 
the peak of a chosen soft X-ray emission line, and then ras- 
tering the instrument pointing to scan an area of interest at a 
fixed wavelength. Line profile measurements were made by 
holding ihe pointing fixed and rotating the crystals to scan a 
region of the soft X-ray spectrum near the wavelength of a 
chosen line. 

The FCS observed the brightest are of the active region of 
interest, at wavelengths around five different emission lines 
Fe xvn (15.01 A), O vm (18.97 A), Mg xi (9.17 A). Ne ix (13 44 
A), and Fe xvm (14.21 A). These lines were chosen because 
they are among the brightest soft X-ray lines within the FC S 
spectral range (about 1.4-22.4 A) that arc produced b> qu ,c! 
active-region plasmas. A peak intensity map and then a me 
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profile observation were made for each line, in the order that 
the lines arc listed. In all. six line profiles were measured, 
because a second Fe xvn line at 15.26 A was included in the 
line* scan near the Fe xvti line at 15.01 A. The observations 
ended with a second Fe xvtt (15.01 A) map. which was only 
partially completed at the end of the SMM orbital day. 

2.1.1. FCS Intensity Maps 

The FCS maps are shown in Figure 1. They range in size 
from 2' to 4’ on a side, with raster step intervals of 10". The 
maps arc of raw counts per second, displayed on a square-root 

scale. , _ . 

Comparison of the first and last maps (both Fe x\n maps) 
shows changes in the active region over the time of FCS obser- 
vation. The region brightened in this approximately 30 minute 
time interval. There were 94 raw counts in the brightest pixel of 
the first Fe xvti map, while the number of counts at the corre- 
sponding position in the last Fe xvn map was about 50% 
hieher. The structure of the active region did not change sig- 
nificantly between the first and last maps. 


2. \ .2. I'CS Speeiral Line Seans 

The FCS line profile observations were made at the position 
of the brightest pixel in the first Fe xvn intensity map. The 
transitions producing each line are listed in Table 1. The O vni 
line results from two different transitions, identified as 
Is 2 5j 2 -lp 2 P Xj2 and \s 2 S { ,-2 p : P 3 The energies of these 
two transitions are so close together that they cannot be 
distinguished from one another in these observations, and so 
they arc treated as one line. The Fe xvnt line is also a result of 
two transitions, identified by Cornille et al. (1992) as 
2p' 2 P X2 -2p A 3d 2 D $ 2 and 2p s 2 P 2 2 -2p A id 2 P } 2 . The wave- 
lengths of the two Fe xvm lines are 14.203 A and 14.215 A 
(Fawcett 1984). To simplify the analysis, these two overlapping 
lines are treated as one line at the average wavelength of 14.21 . 

A line flux is determined by fitting the observed line spec- 
trum with a background level and a Voigtian line profile. The 
assumed line shape is the convolution of the Lorentzian FCS 
instrument response with the Gaussian shape characteristic of 
a Doppler-broadened line. The Fe xvin and O vin lines are 




Ne IX (18:15 UT) 


Fe XVIII (18:21 UT) 


Fe XVII (18:24 UT) 


Fir, I — FCS line peak imcnsiiv maps of AR 4901 between 18:06 and 18:24 IT on 1987 December It. The starting lime of the raster is listed below each map. 
The first Fe xvm 1 15 01 A) raster covers 4' x 4\ The O vim 18.97 A l Mg xi (9 17 Al. Ne ix (13.44 A), and second Fe win 15.01 A rasters com 3’ x 3. and the Fe xvm 
, j 4 •> 1 \ )raMer v * v All the maps are aligned relative to one another and ha\e the same scale, with a step size of 10 " The boxes around each map nave dimensions 
of 4 " v 4 . The last Fe xvn map is incomplete, since the SMM satellite day ended before the raster was completed The images are scaled as the square root of the 
number of raw counts, with a maximum of 100 counts. The number of raw counts in the central pixel of each map is tit first Fe wit map. 94 counts: UlO vm maps. 
34 counts; 1 3) Mg \i map. 28 counts; ( 4 ) Ne ix. 93 counts; (5i Fe win. 23 count*: and 161 second Fe xvn maps. 145 counts. 
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TXBLE 1 

FCS Li^e PkoriiE Observations 


/. Time 1 

Ion I A) (UT) Transition (ergs cm * •' s’ 1 sr~ *) 


Fe xvii 15.01 18:06 :/* ,! S 0 ->' .VP, 5971 ± 238 

Fc xvu 15.26 38.06 Ip" 'S 0 ~2p> 34 'P, 2092 ± 188 

O viii 18.97 38:1? Is *‘5, 2 -2p : P t 2 ; P 3 3 8507 ± 1531 

Mg xi 9.17 18:13 ls : ’Sp-l.Op 233 ± 72 

Nc ix 13.44 18:21 ls : ! S 0 -li 2p , P I 3223 ±183 

Fc win 14.203 18:24 571 ± 257 

14.215 V : P ?2 -2p 4 34 : P 32 


each fit with two Voieiian line profiles, and the measured flux 
is the total flux of the two lines. These line fits are discussed in 
more detail by Waljeski (1992). The line intensities at the Sun 
derived from the FCS line profile measurements are listed in 
Table 1. 


2.2. Broadband Images 

The AS&E High Resolution Imaging Sounding-Rocket 
Payload was launched on a sounding rocket from White Sands 
Missile Range. The AS&E playload was first described by 
Vaiana, Krieger, & Timothy (197?). A more recent description 
is given by Moses. Krieger, & Davis (19S6). The payload con- 
tains a fused-quartz grazing-incidence telescope mirror, five 
interchangeable filters (a white-light filter, two beyllium soft 
X-ray filters, and two aluminized-polypropylene soft X-ray 
filters), and photographic film camera. Full-disk images of the 
Sun on Kodak SO-212 film were obtained for a variety of 
exposure times, ranging from 0.25 s to 60 s. between 18:19:28 
UT and 18:26:38 UT. The entire sequence of observations and 
the properties of each of the filters for this flight are described 
by Waljeski (1992). Two images were selected for this analysis: 
a 3 s exposure through a polypropylene filter, taken starting at 
18:17:08 UT, and a 30 s exposure through a beryllium filter, 
taken starting at 18: 17: 17 UT (see Table 2). 

2.2. 1 . F ilm C alibrat ion 

The film calibration techniques are discussed in detail in the 
doctoral thesis of Waljeski (1992). The film was handled and 
developed as described by Moses et al. (1989). The images were 
digitized with a step size of 2785 x 2285. The film calibration is 
necessarily an iterative process, because the film response is 
wavelength dependent and the incident-filter X-ray spectrum 
depends on the temperature and composition of the observed 
plasma. The shape of the film characteristic curve was found by 
comparing the photographic densities of images with different 
exposure times, as described by Moses et al. (1989). This rela- 
tive calibration technique is based on the histogram cali- 
bration method of Cook, Ewing, & Sutton (1987). The film 
speed (the absolute calibration) is derived from laboratory film 
response measurements at two wavelengths (8 A and 44 A) by 
the iterative method described by Vaiana et al. (1977) and 
Webb & Davis (1985). 


TABLE 2 

AS&E Sol ndi^g- Rocket Payload Observations 



Time 

Exposure Time 

Waveband 

Filter 

(UT) 

(s) 

(A) 

Aluminized polypropylene .. 

38:17:08 

3 

8-39. 44-60 

Beryllium 

18:17:17 

30 

8-29 


2.2.2. Alignment of the Broadband Images with the FCS Observations 

The AS&E broadband images and the FCS observations 
were aligned by comparing the white-light images from each 
instrument to find the position of the solar limb. There are no 
sunspots visible on the SMM white-light image, so the rota- 
tional orientation was fixed by comparing the morphology of 
the broadband images and the first FCS Fe xvn map, which is 
the largest and the brightest of the complete FCS maps from 
this orbit. Figure 2 shows the aligned polypropylene-filtered 
image and Fe xvn map, with the 30" x 30" area of the FCS line 
profile measurements indicated on each image. The accuracy of 
this alignment is estimated to be ± 10" (which is ±3 pixels of 
the broadband image or ± 1 pixel of the FCS map). 

2.2.3. Compensation for the Effects of Scattering 

The point response of the AS&E grazing-incidence telescope 
is sharply peaked with a full width at half-maximum of about 
275. The wings of the point response function are broad, 
because scattering by the telescope surfaces deposits some of 
the incident flux far from the source position. The response of 
the AS&E soft X-ray grazing-incidence telescope and film is 
discussed in detail by Moses et al. (1986). 

A simple scattering correction is made for this data set. as 
described by Waljeski (1992). The effect of the point response of 
the telescope mirror is to scatter the incident energy from the 
coronal sources in the area of observation to points on the film 
outside this area and to scatter energy from coronal sources 
outside the area of observation to inside the area on the film. 
An estimate of the effect of scattering is made by assuming a 
source distribution for a 356" x 356" area and calculating the 
energy scattered into and out of the 30" x 30" area corre- 
sponding to the area of the FCS line profiles. The actual source 
function is not known, but since the point response function of 
the mirror is sharply peaked (if broadly winged), the observed 
image approximates the source. The ratio of the final bright- 
ness in the central area after scattering to the source brightness 
of that area is an estimate of the factor by which the observed 
broadband flux in the area of interest must be chanced to 
compensate for the effects of scattering. 


2.2.4. Selection of the Area of Interest 

To make the intensities measured by the FCS line profile> 
and AS&E broadband images comparable, the broadband 
images are averaged, weighted by the FCS field of view. Thi* 
weighted total broadband flux is then corrected for scattering 
The polypropylene-filtered flux in the region of interest is cor- 
rected for point spread by a multiplicative factor of I S 7 * *' ^ 
the beryllium-filtered flux is corrected by a factor of 1 
(Waljeski 1992). 

The principal sources of uncertainty in the broadband 
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sity measurements arc the film calibration, the alignment with 
the FCS observing area, and the compensation for scattering. 
The total uncertainty in the broadband filtered spectra is esti- 
mated to be about ± 50%. 

3. PLASMA DIAGNOSTICS 

3.1. Temperature and Emission Measure from F CS Lines 
For nearly optically thin emission from an isothermal 
plasma, the intensity emitted at the Sun, /*/.) (ergs cm" s' 
sr' 1 )* in a soft X-ray emission line i at wavelength A, is related 
to the plasma properties of the corona by the expression 

/AA)~PC|(r) n] dl , (1) 

where P is the photon escape probability, t^T) is the emissivity 
of the line (the intensity emitted per unit emission measure at 
temperature T), n t is the electron density, l is the distance 
alone the line of sight, and J n) dl is the emission measure. For 
an escape probability of unity, this expression reduces to the 
optically thin case. 


3.1.1. Emissivity Calculations 

The following expression is used to calculate the emissivity 
of each line: 


[-8.63 x 10-‘<fi>S | 

— A£\ 


L 9 ” Pl 

V kT ), 

n ti n H n. 


( 2 ) 

where A 0 is the wavelength at line center, T is the electron 
temperature. <Q> is the Maxwellian-averaged collision 
strength, B is the branching ratio, g is the statistical weight of 
the lower le^el, A E is the energy difference between the upper 
and lower levels, nJn tX is the ionization fraction, n t fn H is the 
abundance of the emitting element relative to hydrogen, 
is the ratio of the number densities of hydrogen and electrons 
(n a /n r = 0.8 1 , h is Planck’s constant, c is the speed of light, and 
k is Boltzmann's constant. 

The most recent atomic physics calculations available are 
included in the emissivities of the FCS lines. The wavelengths, 
collision strengths, branching ratios, and statistical weights of 
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TABLE 4 

Elemental Ahindances 


Phoiosphcnc Adopied Coronal 

(Grc'cssc & Anders 1991 1 t Meyer 1985b) Coronal 


H 12.00 12.00 12.00 + 0.73 

O 8.93 8.39 9.30 

Ne 8.07 7.54 8.45 

Me 7.58 5.57 8.4g 

f t 7.51* 7.59 8.50 


Mote. — V alues lisied are log of abundance relative to the log abundance of 
H = 1200. 

* The phoiospheric Fe abundance is the more recently determined value 
from the review by Grc'essc et al. 1992. 


the lower levels of each transition are listed in Table 3. The 
single value of QB g listed in Table 3 is the value for the sum of 
the two Fe xvm lines. 

The Maxwellian-averaged collision strength listed in Table 
3. <n>, is the collision strength as a function of energy, inte- 
grated over the electron energies corresponding to a Maxwel- 
lian distribution of velocities. The averaged collision strengths 
are calculated for the velocity distribution of the temperature 
at the maximum emissivity of each line. To be most accurate, 
the Maxwellian-averaged collision strength should be calcu- 
lated as a function of temperature, but this approximation does 
not significantly affect the calculated emissivities in this case. 

The emissivities of the six FCS lines are calculated from 
equation (2) with each of the sets of elemental abundances 
listed in Table 4, the atomic data listed in Table 3, and the 
ionization fractions of Arnaud & Rothenfiug (1985). The 
reasons for choosing these ionization fractions are explained in 
detail in § 5.6. The emissivities for the adopted coronal abun- 
dances of Meyer (1985b) are plotted in Figure 3. 

3.1.2. Optical Depth 

The photon escape probability, P , is included in equation (1) 
to account for optical depth due to resonance scattering. We 



F»C. 3. — Emissivities of '.he lines observed by the SMM FCS, calculated 
from equation (2). with the collision strengths from Table 3 and the ionization 
fractions from Arnaud £i Rothenfiug (1985). The Fe xvn emissivity ploiied is 
that of the line at 15.01 A. The emissivity of the Fe xvn line at 15.26 A has an 
almost identical shape, with slightly lower emissivity values. The broadband 
filtered spectrum integrals for the polypropylene and beryllium fillers arc 
superimposed on the plot for comparison of temperature sensitivity. The fil- 
tered vpectrpm integral has units of ergs cm " : s" l . 


assume that the effect of resonance scattering is only to scatter 
photons out of the line of sight. This assumption is reasonable 
for small optical depths and for the geometry of the observed 
area; since the center of the active region is bright relative to 
the surrounding areas of the corona, the flux scattered out of 
the region of interest should be much greater than the flux 
scattered in. We also assume that the emissivity of the plasma 
is constant along the line of sight in the corona and that the 
line profile shape is Gaussian. In this approximation, the 
photon escape probability is 


P( t 0 ) 


f 1 - CXpj-T 0< ft) 

J J2t<ot 0 


( 3 ) 


where r 0 is the optical depth at line center and <j> is the Gauss- 
ian line profile: 


* “ ex P [- U 2^ ”] ‘ (4) 

The optical depth at line center, t 0 , resulting from resonance 
scattering is given by the expression 


T o 




( 5 ) 


where A 0 is the wavelength of line center and/, u is the absorp- 
tion oscillator strength of the transition from the lower level / 
to the upper level u. The constants e, m, and c are the electron 
charge, the mass of the proton, and the speed of light, respec- 
tively. The line width, <7, is related to the plasma properties by 
the expression 



where M is the atomic mass of the ion and i? w is the non- 
thermal velocity of the plasma. The nonthcrmal velocity is 
assumed to be 20 km s~ l for this analysis. The reasons for 
assuming this value of nonthermal velocity are discussed in 
§ 5.7. The effect of the nonthermal velocity assumption on the 
final results of our analysis are also discussed in § 5.7. 

The intensity ratio of the two Fe xvm lines, which have sig- 
nificantly different oscillator strengths, is a sensitive indicator 
of the optical depth. Combining equation (3) and equation (2) 
shows that the ratio of the fluxes of two lines produced by the 
same ion is independent of abundance and ionization fraction 
and only weakly temperature dependent. For an optically thin 
plasma, the predicted intensity ratio of the two Fe xvn lines 
(the 15.01 A line to the 35.26 A line) is 3.80, the ratio of the 
emissivities. The observed Fe xvn ratio was 1.90 ±0.21. a 


factor of 2 lower than expected for the optically thin c3sc. 
indicating that the photon escape probability is not unity. 

The Fe xvn ratio, calculated as a function of the optical 
depth at line center of the 15.01 A line, is plotted in Figure 4. 
The oscillator strengths used in this calculation./O 5.01 A) = 
2.479 and /(15.26 A) - 0.600, are from Bhatia & Doschck 
(1992). The observed Fe xvn ratio is marked on the cur\c. 
indicating that the optical depth of the 1 5.01 A line is 4.2 ( t 1 
- 1 . 0 ). 

Once the opacity of one line is known, the optical depths o 
all the lines can be calculated. If we assume the Fe xvn ioniza- 
tion fraction and the abundance of iron relative to hxdrogcn. 
the only unknown quantity in the expression for the optua 
depth at line center, given by equation (5). is the co unin 




COMPOSITION OF CORONAL ACTIVE REGION 


9)5 


No. 2. 1994 



density, f n„ dl The optical depth of the 15.01 A line implies the 
f n t d\ - 1.05 x 10 20 cm’ 2 for this region, assuming the 
adopted coronal iron abundance (Meyer 1985b) and J n e dl = 
1.28 x 10 20 cm -2 for the photospheric iron abundance 
(Grevesse et al. 1992). Using these values and the oscillator 
strengths listed in Table 3, the optical depths at line center of 
all the observed lines are calculated and listed in Table 5, for 
both the photospheric abundances and the adopted coronal 
abundances. 

The photon escape probabilities for each of the lines are 
calculated by means of equation (3) and are also listed in Table 
5 For the Fe xvn line at 15.01 A with an optical depth of 4.24, 
the escape probability is 0.358. For this active region, the effect 
of optical depth is significant compared with the uncertainties 
in the intensity measurements for both of the Fe xvu lines, the 
O vm line, and the Ne ix line. . 

Resonance scattering of Fe xvn was also examined for these 
same FCS line data by Schmelz et al. (1992), who found an 
optical depth of 2.4 ( + 0.3, -0.4)for the Fe xvii line at 15.01 A. 
These authors compared the line ratio of Fe xvii (15.01 A)/ 
O vm (18.97 A) to ratios of other FCS lines involving O vm, 
Nc ix. Mg xi, and Fe xvni, under the assumption that only the 
Fe xvii line at 15.01 A is significantly affected by resonance 
scattering. This method underestimates the magnitude of the 
effect if the O vm line is also significantly affected, as shown 
here. The present study, w-hich compares two affected Fe xvn 
lines to obtain a more direct, abundance-independent measure 
of the resonance scattering, provides a more self-consistent 


treatment of the scattering opacity of all the lines, both in the 
calculation ol optical depths and in the determination of 
plasma parameters. 

3.2. Temperature and Emission \f easurefrom 
Broadband Images 

For an isothermal optically thin plasma, the intensity F i% in 
ergs cm -2 s’ 1 , incident on the film at the focal plane of the 
broadband telescope, through a given filter, i is related to the 
temperature, emission measure, and composition of the emit- 
ting coronal plasma by the equation 

F, = J fUK J ", J dl , (7) 

where A is the geometrical collecting area of the telescope, /is 
the telescope focal length, is the product of the filter trans- 
mission and the telescope reflectivity as a function of wave- 
length, and p(A, T) is the power emitted at wavelength /. at 
temperature 7\ per unit emission measure. 

3.2.1. The Broadband Filtered Spectra 

For an isothermal plasma, the integral over wavelength in 
equation (7) (the filtered spectrum integral), expresses the tem- 
perature dependence of the flux equation. The calculation of 
the filtered spectrum integrals is described in detail by Waljeski 
(199**) 

The solar spectrum, p( A, T), is calculated between 8 and 60 A 
usine a program written by Lemen (1990). In this wavelength 
range, both continuum emission and a large number of emis- 
sion lines (about 500) produced by many different elements and 
ionization species contribute to the solar spectrum. Lemen s 
program generates a model solar spectrum using the line emis- 
sivities of Mewe, Gronenschild, & van den Oord (1985) and a 
model of the solar soft X-ray continuum. The emissivities tabu- 
lated by Mew'e et al. (1985) are not the most recent calculations 
for many of the lines, but they are in a form that is convenient 
for the broadband data analysis and are sufficiently accurate 
for this purpose, considering the large number of lines that 

contribute to each of the broadband filtered spectra. 

The relative contributions of lines produced by all ions of 
Ne, O, Mg, and Fe to the filtered solar spectra are plotted in 
Ficures 5a and 5b for the adopted coronal abundances of 
Meyer (1985b). The percentage contribution of the continuum 
emission is also plotted. This plot shows that the filtered 
spectra are dominated by line emission, particularly from i f0n 
ions, at active-region temperatures^ about 2x10 K to 
3 x 10 6 K. The percentage contribution of iron lines to the 
filtered spectra is about 80%, and the percentage contribution 
of the continuum is less than 20%, over most of this tem- 
perature rage. 


TABLE 5 

Optical Depths and Photon Escape Probabilities 


Photospheric Abundances Adopted Coronal Abundances 


Ion 

A o 

(A) 

*0 

P 

*0 

P 

Fc xvu 

Fe xvn 

O vm 

Mg xi 

15 01 
15.26 
18.97 
9.17 

4.24 

1.04 

7.59 

0.80 

0.358 

0.716 

0.233 

0.770 

4.24 

1.04 

1.82 

0.65 

0.358 

0.7)6 

0.581 

0.806 

0.797 

0.932 

Ne ix 

Fe xviii 

13.44 

1421 

2.79 

0.20 

0.466 

0.932 

0.69 

0.20 
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Fig. 5. — Percentage contribution of all the lines of Fe. O. and Ne and of the continuum, to the filtered solar spectra. The contributions of lines of Me to the 
filtered spectra are less than 1 % for most temperatures, and so are not plotted, (a) Percentage contributions of each element to the aluminized-polypropylene-fihered 
soft X-ray spectrum (8-39. >*4-60 A); and (M Percentage contributions to the beryllium-filtered soft X-ray spectrum (S--9 A). For both filters, the Fe lines dominate 
the filtered spectrum over a large range of temperatures, making up about 80% of the lota) flux through each filter at temperatures characteristic of quiet coronal 
active regions (about 3 x 10 6 K). The total contribution of all the elements is not equal to 100% at all temperatures because other elements, including C. C a, Ni. N. 
and S. also contribute to the filtered solar spectra in this temperature range. 


The broadband spectrum is assumed to be optically thin. 
This assumption is made because it is impractical to calculate 
the optical depth and photon escape probability for each of the 
more than 500 lines that contribute to the filtered spectra, and 
because this assumption is not likely to have a significant effect 
on the integrated filtered flux (relative to the 50% uncertainty 
in the broadband intensity measurement). Only the brightest 
lines are likely to have large optical depths in this case, so most 
of the lines that contribute to the broadband spectrum should 
not be significantly affected by resonance scattering. 

The filtered spectrum integrals for the polypropylene and 
beryllium filters are superimposed on the line emissivity plot 
(Figure 3) for comparison of the temperature responses. The 
units of the filtered spectrum integrals (ergs cm' 2 s' 1 ) are not 
the same as the units of the line emissivities (ergs cm' 3 s' 1 
sr' 1 ) because the waveband is selected by the telescope reflec- 
tivity and the filter responses, so the broadband instrumental 
response can not be separated from the emissivity of the pass- 
band (Waljwski 1992). 


4. RESULTS 

The relative abundances of iron, oxygen, magnesium, and 
neon in the corona are determined by finding a model active- 
region plasma with a temperature, emission measure, and ele- 
mental composition that would produce the observed line and 
broadband intensities. We then use the additional information 
provided by the optical depths of the Fe xv» lines and the 
geometry of the region as a diagnostic of absolute abundance. 

4.1. Relative Abundances of Fe, O, Mg, and Ne 

We first consider the simplest temperature and emission 
measure distribution: the isothermal model. If the observed 
active region is composed of isothermal plasma, then it should 
be described by a single temperature and emission measure. 
The emission measure. EM (T) = J n 2 dl , necessary for a plasma 
of a given temperature and composition to produce the 
observed flux is found by inverting the intensity equations (for 
FCS lines, eq. [1], and for broadband intensities, eq. [7]). The 


derived emission measures are plotted as a function of the 
assumed isothermal temperature in Figure 6 a for the photo- 
spheric abundances of Grevesse & Anders (1991) and in Figure 
7a for the adopted coronal abundances of Meyer (1985b). 
These plots provide a qualitative indication of how well the 
isothermal model plasma fit the observed intensities. 

A quantitative measure of the consistency of the observed 
intensities with the plasma model is provided by the reduced y 2 
of the expected intensities for each temperature and emission 
measure: 


r(7\ EM) 


1 f f /,- - 1JJ\ EM) ~P 
•V - 1 L e, 


(?) 


where N — 1 is the number of degrees of freedom (i.e., the 
number of observations minus one), f is the observed intensity. 
Ij(T 9 EM) is the intensity predicted from an isothermal model 
plasma with temperature T and emission measure EM, and c, 
is the standard deviation of each observed intensity. 

The temperatures of isothermal plasmas consistent with the 
observations are found by the following method: The reduced 
y 2 (7\ EM) is calculated for model isothermal plasmas with a 
range of temperatures (1 x 10 6 K to 10 x 10 6 K). At each 
temperature, we find the model plasma that best reproduces 
the observed intensities, by identifying the emission measure 
that yields the minimum value of y 2 (T, EM) at that tem- 
perature. A minimum reduced /o(7~) of unity or less indicates 
that an isothermal plasma with the given temperature is con- 
sistent with the observations for some emission measure or 
range of emission measures. 


4.1.1. Photospheric Abundances 

In Figure 6a , we have plotted the emission measures derived 
for each line and broadband intensity measurement as a func- 
tion of the assumed temperature for a model plasma having 
photospheric abundances (Table 4). On this plot, the emission 
measures determined from all the line and broadband observa- 
tions do not appear to be mutually consistent for an isotherm-* 
plasma of any temperature. The FCS lines from two clusters on 
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Fig. 6 a 


Fig. 6 b 


riw. , . 

emitting region. Each Curve represents the intensities and from equation (7) for the broadband intensities, (b) Minimum (with respect to 

emission measure is calculated from equation (1) for isothermal plasma with photospheric abundances, as a function of temperature. The minimum 

»"S 5 «ta entire temperature range, indicating .ha, the observed line and broadband intensities are not consistent wuh the 

emission from an isothermal photospheric plasma. 


the temperature-emission measure plot: one cluster near the 
emission measure determined by the broadband 
polypropylene-filtered spectrum, and another cluster at an 
emission measure about an order of magnitude higher. The 
croup with lower emission measure includes the lines of 
Fe xvii, Fe xvm, and Mg xi, all low-FIP elements, as well as 
the broadband filtered spectra. As shown in Figures 5a and 56, 
the broadband filtered spectra are dominated by emission lines 
of iron, which is a low-FIP element. The higher emission 
measure eroup contain the lines of O vii and Nc ix, both high- 
FIp elements. The separation of the emission measures deter- 
mined from the lines of high-FIP and low-FIP elements for the 
photospheric abundances suggests the fractionation of the 


composition of the solar corona with FIP discussed by Meyer 
(1985b). , 

The minimum reduced x of the intensities, plotted in Figure 
66, is not less than 8.5 for any temperature. Therefore, the 
probability that the joint line and broadband observations are 
consistent with an isothermal plasma of photospheric com- 
position for any temperature or emission measure is less than 
0.1% (Bevington 1969). 

4.1.2. Adopted Coronal Abundances 

The emission measures derived from all of the observed 
intensities are plotted in Figure 7 a for an isothermal plasma 
with the FlP-fractionated adopted coronal abundances of 




Fig. 7. — (j) Emission measure as a function of temperature, assuming the FlP-biased temperature, 

measure distribution. Each curve represents the emission measure derived from ^ measured 1 line or inlensilies . {b) Minimum reduced ? of the 

The emission measure is calculated from equation (!) for the FCS *ne i n< j abundances of (Mcver l9S5b). as a function of temperature. The value of 

WSSSrKSS 5 v - — in ihis — ** 

l he emission from an isothermal plasma having the adopted coronal relative elemental abundances 
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Mc\er i ] 9Srb ». On this plot, the emission measures determined 
by the joint observations seem to be consistent with one 
another for activ e*recion temperatures (between 2 x 10 6 K 
and 3 x 10* K l 

The minimum reduced chi-squared of the intensities, plotted 
in Figure 7fi. is less than one over a range of temperatures from 
2.5 x 10* K. to 3.0 x 10 6 K, indicating that these data are 
consistent with an isothermal plasma of adopted coronal com- 
position, having any temperature in this range. The ranee of 
values of emission measure corresponding to this temperature 
range were 3.0 x ]0 Z9 cm' 5 to 1.7 x 10 29 cm" 5 . The 
minimum value of /1{T) corresponds to a probability of fit of 
90% (Bevington 1969) for an isothermal plasma with a 
temperature of 2.8 x 10 6 K and an emission measure of 
2.3 x 3 0 ;s> cm" 5 . 


4.2. Absolute Abundances 

We have confirmed, assuming an isothermal plasma, that 
the ratio of the abundances of the high-FlP elements (O and 
Ne) to the low-FIP elements (Mg an 6 Fe) is about a factor of 
3.5 less in the corona than in the photosphere. The additional 
information provided by the optical depths of the Fe xvn lines 
and the geometry of the region makes it possible to constrain 
the abundances of the heavier elements relative to the abun- 
dance of hydrogen, as well as to determine the electron density. 
Assuming that the plasma is uniformly distributed along the 
line of sight through the active region, the active-region density 
and path length can be determined from the emission measure, 
j n 2 dl and the line-of-sight integral of density, J n 4 dl , using the 
relations 


and 




,f "id 



(9) 


(10) 


The abundance of iron relative to hydrogen may then be deter- 
mined by comparing the path length calculated from the 
optical depth of the Fe xvn line to the geometry of the active 
region seen in the broadband images. Examining equation (5), 
we see that the expected path length is inversely proportional 
to the absolute abundance of iron, n T Jn H . We assume that all 
the parameters in the equation, including the ionization equi- 
librium and the atomic parameters are known. Therefore, it is 
possible to determine the absolute abundance of iron, n T Jn H , 
given the measured optical depth and the path length. 


4.2.1. Density and Path Length for Meyer Abundances 

The active-region density and path length are first calculated 
for the adopted coronal abundances of Meyer (1985b). The 
value of the integral of density along the line of sight derived 
from the ratio of the two Fe xvn lines, J n t dl = 1.05 x 10 20 cm, 
and the emission measure derived from the observed inten- 
sities, jn*dl = 23 x 10 29 cm" 2 ( + 0.8 x 10 29 cm 2 , 
— 0.4 x 10 29 cm" : ) give a density, n t =* 2.2 x 10 9 cm" 3 
( + 0.8 x 10 9 cm" 3 , -6.4 x 10 9 cm" 3 ). This density is in rea- 
sonable agreement with densities derived from density- 
sensitive lines ofFe xii, Fe xm, Fe xiv, and Fe xv(Dere 1982). 

The path length obtained for these abundances is 4.7 x 10 10 
cm (-1.8, +0.9) or 653" (-281", +228"). This path length is 
much longer than could be expected for this region. The broad- 


band image in Figure 2 shows that the emission from the entire 
active region spans only about 150”, while the calculated path 
length is about a third of the diameter of the Sun. It is difficult 
to sec how such a long path length could be associated with 
this active region. This active region is near the limb of the Sun 
(at about 70‘ west), but the 30" x 30" area of interest is entirely 
on the disk in this line of sight, so material from the other side 
of the limb could not be extending the path length. 

The path length along the line of sight cannot be measured 
directly from the broadband or the FCS images, but a reason- 
able estimate may be made by measuring the characteristic 
lengths of the region. The largest dimension (a southeast-to- 
northwest diagonal) of the entire active region is about 150". 
The characteristic path length is estimated from profiles 
through the area of interest, made approximately east to west 
across the maps of Jn 2 J/ and (\n]dl)' 2 constructed from the 
broadband images. From the measured full width at half- 
maximum of these profiles, we estimate that the path length is 
80" ± 40". 

4.2.2. Determination of Absolute Abundances 

The path length of 653" derived using the adopted coronal 
abundance of Fe (Meyer 1985b) is clearly inconsistent with the 
observed dimensions of this active region. The derived and 
observed path lengths can be made consistent if we take hydro- 
gen to behave as a high-FIP element, as suggested recently bv 
Meyer (1993a, 1993b). In order to account for the observed 
path length, the abundance of iron relative to hvdrosen must 
be 3.20 x 10‘ 4 (-2.12 x 10' 4 , +2.00 x 10 _i ). This iron 
abundance is higher than the value given by Meyer (1985b) b> 
a factor of 8.2 ( — 5.4, +5.0). 

The calculated temperature, electron density, optical depth, 
and escape probability are not affected by this change in the 
assumed value of the hydrogen abundance, while the emission 
measure, and the column density are inversely proportional to 
n T Jn H . For the absolute abundances determined here. 

J n 4 dl = 2.8 x 10 28 and J n t dl = 1.3 x 10 19 . 

Since the relative abundances of Fe, O, Mg. and Ne giv en bv 
Meyer (1985b) are consistent with these observations, wc are 
able to calculate the abundances of O, Mg, and Ne relative to 
hydrogen, given the calculated value of n Fe /n H . These absoluic 
abundance values are listed in the last column of Table 4. 

If we assume that all the relative abundances of the heavier 
elements given by Meyer (1985b) for the solar corona arc 
correct, a complete set of coronal abundances that is consistcni 
with the present results can be obtained. These new coronal 
abundances, listed in Table 6, are calculated by increasing the 
Meyer (1985b) adopted coronal abundances of all the elements 
except hydrogen by a factor of 8.2. 

The uncertainty in the hydrogen abundance listed in Tables 
5 and 6 results from the measurement of the Fe xvn intensities 
and the path length. The simplicity of our analytical models, 
which assume homogeneous conditions along the line of sight, 
may contribute systematic errors, in addition to the stated 
measurement errors, if the material along the line of sight is not 
uniformly distributed. 

4.2.3. The Effect of Inhomogeneous Density 

To calculate the absolute abundances, we assumed a 
uniform density distribution in deriving the path length. A*, 
from equations (9) and (10). In this section, we consider the 
effect of assuming an inhomogeneous density distribution on 
the expected path lengths and abundances. 
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TABLE 6 

Coronal Abundances 


Element Abundance 


H 12.00 ±0.73 

He 11.90 

C 9.28 

N 8.50 

O 9.30 

Ne 8.45 

Na 7.35 

Mg 8.48 

Al 7.35 

Si 8.50 

S 7.84 

Ar 7.24 

Ca 7.38 

Fe 8.50 

Ni 7.24 


Note. —Values listed are log of abundance normalized 
so that the log abundance of H ** 12.00. The university 
listed is the uncertainty in the abundance of the heavier 
elements relative to the abundance of H. 

Two reasonable models of nonuniform density are con- 
sidered: a hydrostatic model and a filamentary model (Fig. 8). 
These models are functions of a characteristic length h and 
density n 0 . The hydrostatic distribution reproduces the decline 
in intensity with height observed in the broadband X-ray 
image. In the hydrostatic model, h is the scale height of the 
distribution. The filamentary model has a filling factor of 
and approximates the looplike structure of active regions. 




Distance (1) 

Fig. 8. — Two inhomogeneous density distributions characteristic of active 
regions. These models are functions of a characteristic length h and density n 0 
and are scaled so that the emission measure is equal to nj 8. 
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Each density distribution is scaled so that the emission 
measure is ec ual to h. 

For the hydrostatic model, the path length Af calculated 
from equations (9) and (10) equal to 2h. Consequently, if we 
had assumed a uniform density distribution with a character- 
istic length h y we would have underestimated the path length 
by a factor of 2 and thus overestimated the abundance of iron 
relative to hydrogen by a factor of 2. For the filamentary 
model, the effect of nonuniform density on the absolute abun- 
dances works in the opposite sense : For a filling factor of the 
calculated path length is h/ 2, so if a uniform density were used, 
the path length would be overestimated and the iron abun- 
dance underestimated by a factor of 2. The actual density dis- 
tribution in this active region is probably some combination of 
a gradually decreasing density with looplike structures, and the 
effects of the two density distributions on the derived abun- 
dances w'ould tend to cancel. 

4.3. Interpretation of the Results 

We have demonstrated, for this active region, that the 
low-FIP elements are overabundant in the corona by a factor 
of 8.2 (-5.4, +5.0) relative to hydrogen, while the high-FIP 
elements other than H arc overabundant in the corona by a 
factor of 2.3 (—5.4, +5.0). The uncertainty in our measured 
hydrogen abundance limits a determination of whether the 
abundances of the high-FIP elements (except H) are nearly the 
same in the photosphere and corona or are enhanced. Equiva- 
lently, this uncertainty prevents us from establishing whether 
hydrogen fractionates in the corona to the same degree as the 
other high-FIP elements. However, it is clear, in spite of the 
magnitude of this uncertainty, that hydrogen does not behave 
as a low-FIP element in the corona. 

5. DISCUSSION OF ANALYSIS TECHNIQUE 

Next, we consider to what extent our relative and absolute 
abundance measurements might be affected by various factors, 
including the instrumental calibration, the evolution of the 
region over time, and the effects of optical depth, a multi- 
thermal plasma, the assumed iron ionization balance, or the 
assumed line broadening. We show that these factors have no 
significant effect on our results. Specifically, these additional 
considerations do not change the measured relative elemental 
abundances cr the conclusion that hydrogen does not fraction- 
ate as a low-FIP element in the corona. 

5.1. Instrumental Calibration 

If the relative instrumental responses of the AS&E 
sounding-rocket payload and the FCS were not correctly 
understood, vc would expect to find two separate solutions of 
emission measure (for the same isothermal model) from the 
observations, one from the FCS line data and one from the 
broadband images. 

For an isothermal plasma with Meyer’s adopted coronal 
abundances, ihe emission measures from the broadband data 
are about 30% low'er than the emission measures derived from 
the FCS lines at active-region temperatures. This discrepancy 
is within the uncertainties of the data, as indicated by the 
minimum reduced * 2 values of less than unity in this range of 
temperatures. 

For an isothermal plasma with photospheric composition, 
the FCS line intensities are inconsistent with one another as 
well as with broadband observations, indicating that problems 
with the relative instrumental calibration could not have the 
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effect of creating the apparent FIP bias in the coronal abun- 
dances. 

5.2. Evolution of the Active Region 

As discussed earlier, active region AR 490) was observed to 
have brightened in the Fe xvn line at 15.01 A: the central pixel 
of the last Fe xvn map, the area of interest, had about 50% 
more counts than the corresponding pixel of the first Fe xvn 
map. These two Fe xvn maps span the time of FCS observa- 
tion, so this brightening could have happened at any time 
during this approximately 50 minute period. It seems unlikely 
that an increase in the intensity of the region caused the appar- 
ent FIP bias in the coronal abundances. The O vni and Ne ix 
measurements were taken near the middle of the observing 
period, at 18:13 UT and 18:21 UT, respectively. These two 
lines of high-FIP elements are not consistent with the lines of 
low-FIP elements, for photospheric abundances. A change in 
the brightness of the region of the magnitude of the FIP differ- 
ential (a factor of 3.5) would be apparent in the Mg xi data, 
which were taken simultaneously with the O vm data, as well 
as in the broadband images (at 18:17 UT) taken in the time 
period of the O vm (18:08-18: 15 UT) and Ne ix (18:15-18:21 
UT) observations. Since such a decrease is not observed in the 
Mg xi fiux or in the broadband fiuxes, the FlP-differential of 
the elemental abundances deduced from these data is not likely 
to be an effect of a change in the brightness of the active region 
over time. 

5.3. Optical Depth 

The optical depths of some of the lines, especially the Fe xvn 
line at 15.01 A, are significant compared with the uncertainties 
of the flux measurements, so the effects of optical depth must 
be taken into consideration in the determination ofihe coronal 
abundances. If we assumed that these FCS lines were optically 
thin, then the observed line and broadband intensities would 
be inconsistent with a plasma having either set of relative 
abundances, at any temperature: the minimum reduced y 2 of 
the intensities of an optically thin plasma is 6.2 for the adopted 
corona] abundances (Meyer 1985b) and is 11.6 for photo- 
spheric abundances. 

The optical depth of the lines gives us the additional infor- 
mation needed for a diagnostic of absolute abundance. It is not 
possible from our temperature and emission measure models 
alone to determine whether the abundances of the high-FIP 
elements are depleted in the corona relative to the photo- 
sphere, as proposed by Meyer (1985b), or whether the abun- 
dances of the low-FIP elements are enhanced, as recent 
evidence suggests (Meyer 3993a, 1993b), because there is no 
separate measurement of the continuum in this data set from 
which to directly measure the abundance of hydrogen. For the 
adopted coronal abundances (Meyer 1985b), emission from 
iron lines dominates the broadband spectra at active-region 
temperatures, as show n by Figures 5 a and 5b. The continuum 
contributes less than 20% of the total filtered spectrum of both 
of the broadband filters. Thus, if all the lines in this region 
actually were optically thin, then no information on absolute 
abundances would be present in the data. 

5.4. Broadband Temperature Map 

The temperature distribution within the 30" x 30" area of 
interest can be examined by looking at a broadband tem- 
perature map made from the ratio of the images through two 
different filters. Figure 9 shows a contour plot of the tem- 



Fjg. 9.— Temperature msp. in uniis of 10* K, of lhe area of joint observa- 
tion near ihc center of AR 4901. The temperature map is calculated from the 
ratio of the polypropylene-filtered and the beryllium -filtered broadband soft 
X-ray images, assuming that the region is isothermal along the line of sight. 
The adopted coronal abundances (Meyer 1965b) are used in the calculation. 
The line-of-sight temperatures in this region range from 2.3 x 10* K. to 
3.0 x 10 6 K. 


perature derived from the broadband images in the area 
observed by the FCS collimator. (This is the area over which 
the line profile observations wxre made.) Within the area of 
interest, the line-of-sight temperatures (for 6" x 6" resolution 
elements) vary from 2.3 x 10* to 3.0 x JO 6 K. This range of 
temperatures corresponds to values of the minimum reduced 
y 2 of near unity for the Meyer (1985b) adopted coronal 
abundances (Figure 7 b). Most of the line-of-sight temperatures 
in this map are within the range of isothermal temperature 
solutions (2.5 x 30 6 to 3.0 x 10 6 K), especially considering 
that the pyramidal field of view of the FCS collimator weights 
the material at the center more heavily in the intensity mea- 
surement. 

5.5. Muhithermal Analysis 

Next we consider the effect of a multitemperature plasma on 
the abundance diagnostic. We will allow r the active-region tem- 
perature to vary along the line of sight as w ell as perpendicular 
to it. For a nonisothermal plasma, the definition of emission 
measure, $n 2 dl, is no longer applicable. In the multithcrmal 
case, we use the differential emission measure, ;{T). 

Empirically determined differential emission measures 
(OEMs) (e.g., Dere & Mason 3993) can be w’ell approximated 
at corona] temperatures by a Gaussian shape. We assume that 
the DEM for this active region, C(T), has a Gaussian shape 
described by three parameters — the temperature correspond- 
ing to the peak of the Gaussian curve ihe half-w idth of 

the curve (AT), and a multiplicative factor (f 0 ): 

_ - f (r-r m „) 2 ] ll!f 

«r)~t.«p[- :iif — J- 

For a differential emission measure of this form, the intends ^ 
a coronal emission line at the Sun is given by the exprew ion 
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Similarly, the broadband flux incident on the film is given in 
terms of the differential emission measure by the expression 


F,= 


A 


4it/- 



(T - T m J : 
2(AT) : 


p(/., log T)>i/i/.)d/.dT . 


(13) 


The inteeral over temperature in equations (12) and (13) is 
evaluated over the entire range of formation temperatures for 
each line or broadband spectrum. 

The value of the minimum reduced jr is calculated^ for a 
ranee of values of the peak temperatures T m „ (1 x 10 K to 
10 x 10 6 K) and half-widths, AT (0.1 x 10 K to 10 x 10^ K). 
The scaling factor c 0 that yields the minimum reduced x tor 
each (7 ,~AT) pair is determined. If Zo(F) is less than or equal 

to unity for any value of the DEM scaling factors, then the 
data are consistent with a plasma having a Gaussian DEM 
described by the given T miX and AT. ncu 

The values of the minimum reduced y- for each pair of DhM 
parameters, T mis and AT, are plotted on a contour plot in 
Figure 10a for a plasma with photospheric abundances and in 
Figure 10b for one with the adopted coronal abundances. As in 
the isothermal case, there is no Gaussian DEM model for 
which yl(T) is equal to or less than eight, indicating that the 
observations are not consistent with a plasma having photo- 
spheric abundances and a Gaussian DEM. For the adopted 
coronal abundances (Meyer 1985b), shown in Figure 10b, there 
are a range of OEMs for which the *g(T) is less than or equal 
to unity. This indicates that the data are consistent with multi- 
thermal plasmas having the adopted coronal relative abun- 
dances (Meyer 1985b) and a range of Gaussian DEMs. 

The dotted lines plotted over the contours in Figure 10b 
show the isothermal plasma temperatures that are consistent 
with the data. For narrower Gaussian DEMs, the range of 
peak temperatures that are consistent with the data 
approaches the range of isothermal temperatures solutions. 
Plasmas having Gaussian DEMs with half-widths of up to 
6 x 10 5 K are consistent with the data for peak temperatures 
near to or lower than the temperatures of the isothermal solu- 
tions. It is not possible to distinguish between the multithermal 
models and the isothermal models on the basis of these data 


any plasma model that fits the data satisfactorily (i.e., has a 
minimum reduced y 2 of one or less) is a possible solution for 
the plasma properties of the region. 

We have shown that there are multithermal DEM models as 
well as isothermal models that arc consistent with the observa- 
tions for plasmas with Meyer’s (1985b) adopted coronal rela- 
tive abundances, but not for those with photospheric 
abundances. Therefore, it is not likely that the apparent deple- 
tion of the high-FIP abundances relative to the low-FIP abun- 
dances in the corona is an effect of a multithermal temperature 
distribution. 


5.6. Iron Ionization Fractions 

The ionization balance calculations used in this paper are 
those of Amaud & Rothenflug (1985). Recently, the ionization 
equilibrium of iron has been recalled by Arnaud & Raymond 
(1992). Saba & Strong (1992) have found difficulties in the 
interpretation of FCS data with these new iron ionization frac- 
tions. Since we have measurements of only three iron lines (two 
Fe xvit lines and one Fe xvm line), it is not possible, given the 
uncertainties in the measured fluxes of these lines, to determine 
which of these iron ionization balance calculations is more 
accurate. In our analysis, we chose to use the Arnaud & 
Rothenflug (1985) ionization fractions for consistency, as they 
are assumed in the emissivities (Mewe et al. 1985) used in the 
interpretation of the broadband spectra. In this section we 
discuss how the assumption of the Arnaud & Raymond (1992) 
iron ionization fractions w’ould affect our results. 

5.6. 1 . Relative Abundances 

The relative abundance results are the same for the assump- 
tion of either iron ionization balance — the observations are 
not consistent with an isothermal plasma of photospheric 
abundances, but are consistent with an isothermal plasma with 
the adopted coronal abundances of Meyer (1985b) for a range 
of temperatures from 2.6 x 10 6 K to 3.0 x 10 6 K. 

5.6.2. Absolute Abundances 

The .absolute abundance calculation is affected by the choice 
of iron ionization fractions, but the result that hydrogen does 
not behave as a low-FIP element in the corona is not changed. 
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Fic. 10a 

Fig. 10. — Contour maps of minimum reduced showing the fit of Gaussian 
abundances and for ibt ihe adopted coronal abundances (Me>er 1985b). 



Fig. 10b 

DEM models to ihe data as a function of T m ,, and ST for (a) the photospheric 
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For ihc Arnaud & Raymond (1992) values, the path length is 
about 50% longer than the path length calculated in § 4.2. This 
leads to an iron abundance relative to hydrogen of 7.4 x 10" 4 
(~4.9 x 10 4* 4.6 x 10'*), meaning that the low-FIP ele- 

ments must be overabundant in the corona by a factor between 
6 and 31. 

5.7. Somhermal Velocities 

The statistics of these FCS line profile data do not allow' us 
to accurately measure the excess line broadening. Therefore, 
for this analysis, a value of the nonihermal velocity was 
assumed based on other coronal observations. Measurements 
of nonihermal broadening of coronal active-region lines vary 
from 10-30 km s' 1 (Cheng, Doschek, & Feldman 1979) from 
cooler EUV lines, to 26 ± 5 km s' 1 from a 1.3 x 10 6 K Me 
line (Hassler et al. 1990), to44kms _1 for the average of a 
3 x 10 6 K active region (with a range 0-60 km s' 1 over the 
active region) from the Mg xj line at 9.17 A (Saba & Strong 
1991). 

We arbitrarily assumed 20 km s" 1 as an intermediate value 
of nonihermal velocity and then investigated the effects of 
assuming nonihermal velocities ranging from Okms" 1 to 60 
km s“ 1 on our abundance results. The nonthermal velocities of 
all the FCS lines were assumed to be same, although only the 
nonihermal velocity assumed for the Fe xvn line at 15.01 A has 
a significant effect on our results. 

The relative abundances determined from these observa- 
tions are unchanged for the assumption of any nonthermal 
velocity in this range. 

The effect of the nonthermal velocity assumption on the 
absolute abundance of iron relative to hydrogen derived from 
observations of this active region is dependent on the magni- 
tude of the excess velocity relative to the Fe xvn (15.01 A) 
thermal velocity of 29 km s" 1 : if w>e assumed that the lines 
were thermally broadened only (as Schmelz et al. 1992 did for 
these same FCS observations), then the derived absolute abun- 
dance of iron would decrease by only 20% from the value 
obtained in §4.2.2 for r nt «= 20 km s* 1 . For a nonthermal 
velocity of 29 kra $~\ the derived absolute abundance of iron 
is 30% greater than it is for v nt = 20 km s' 1 . However, the 
absolute abundance of iron determined from these observa- 
tions is significantly increased over our previous result of 
§4.2.2 if nonthermal velocities larger than 29 km s” 1 are 
assumed : for r m *= 40 km s ’ 1 the absolute abundance of iron is 
increased by a factor of 2, and for r nt « 60 km s“ 1 the absolute 
abundance of iron is increased by a factor of 3.7. Qualitatively, 
the effect of the assumed nonthermal broadening on this 
analysis does not change the conclusion that the abundance of 
iron is enhanced in the corona relative to hydrogen and that, 
therefore, hydrogen behaves more as a high-FIP element than 
as a low’-FIP element in the corona. Qualitatively, the degree 
of iron enhancement in the corona is significantly greater for 
the assumption of nonthermal contributions to the line widths 
of greater than 30 km s' L 

6. SUMMARY 

The simultaneous broadband and line intensity data from 
the center of the observed coronal active region confirm the 
FIP separation of Meyer’s (1985b) adopted coronal abun- 


dances of O, Ne, Mg, and Fe and put these abundances on an 
absolute scale with respect to hydrogen. The low-FIP heavy 
elements are found to be overabundant in this coronal active 
region by a factor of 3.5 relative to the high-FIP elements. We 
find that hydrogen behaves more like a high-FIP element in 
the corona, as suggested by Meyer (1993a, 1993b) and 
Feldman (1992), than like a low-FIP element (Meyer 3985b): in 
this coronal active region, hydrogen is a factor of 8.2 (-5.4, 
+ 5.0) less abundant relative to the heavy dements than in the 
photosphere. 

The systematic uncertainties in this analysis arise primarily 
from the assumptions of the ionization fraction, the density 
structure of the active region, and the nonthermal line 
broadening. These uncertainties do not change our relative 
abundance results, but they all tend to increase our derived 
absolute abundance of iron. Considering the systematic and 
statistical uncertainties we conclude that the abundance of 
hydrogen relative to iron must be depleted in the corona com- 
pared to the photosphere by a factor of at least 2.3 in this activ e 
region, but it may be depleted by more than an order of magni- 
tude more. Even in the most extreme case, the systematic 
uncertainties do not affect the conclusions that, for this active 
region, hydrogen does not fractionate as a low-FIP element in 
the corona and the relative abundances of O, Ne, Me, and Fe 
conform to the FlP-separation discussed in the review bv 
Meyer (1985b). 
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